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TEMPERATURES below approximately 10°K can _ be 
obtained only through the use of helium gas as a refriger- 
ant. Moreover, in order to liquefy helium it is necessary to 
compress the gas to approximately 200-300 Ib/in?, cool it 
to approximately 15°K or colder, and expand it through 
a Joule-Thomson cycle to | atm or less. All liquefiers or 
refrigerators must necessarily have this feature in common. 
Refrigerators differ from each other according to the way 
in which the initial cooling from 300°K to 15°K is 
accomplished. 

The classical method, and the one first used to liquefy 
helium over fifty years ago, utilized the boiling of liquefied 
nitrogen and liquefied hydrogen under reduced pressure as 
the method of precooling. For liquefying batches of helium 
this method is still used in some laboratories which have 
an adequate supply of liquid hydrogen. However, it is not 
an attractive method for producing refrigeration con- 
tinuously in a closed cycle. z 

A considerably improved method ot pioviding precool- 
ing without using any refrigerants other than helium was 
developed by Collins,” and has led to the widely accepted 
ADL-Collins laboratory liquefier.* In this system, portions 
of compressed helium gas are expanded in expansion 
engines, following Kapitza,* at intermediate temperatures 
to produce the refrigeration required for precooling the 
Joule-Thomson stream. The Collins method represents a 
great step towards the realization of a closed cycle refrig- 
eration system because of its simplicity, safety, and in- 
herent efficiency. However, there are certain problems of a 
practical nature inherent in the engine design which tend 
to limit the mechanical reliability. Foremost among these 
is the requirement for a closely fitting piston—cylinder 
combination (which, of course, cannot be lubricated), 
and the necessity of having mechanically operated engine 
valves which operate at very low temperatures. 
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More recently, adiabatic expansion machines of a quite 
different sort have been developed which could conceivably 
serve as precooling devices for helium liquefiers. The early 
concepts of Kirk and Stirling in which a reciprocating-flow 
thermal regenerator is used in place of a countercurrent 
heat exchanger have recently been refined and developed 
effectively into an air liquefier by the Philips Company.® 
The system is characterized by the use of two articulated 
pistons which provide a compression space and an expan- 
sion space, interconnected through a thermal regenerator. 
Unlike the Kapitza or Collins system, there are no valves 
required either for the compression or expansion sections 
of the Philips apparatus, and this feature may be expected 
to be advantageous from the point of view of mechanical 
reliability. On the other hand, the inability to separate 
physically the expansion and compression sections of the 
machine imposes a considerable degree of inflexibility and 
almost eliminates the possibility of using such a refriger- 
ator as a multi-stage precooler. 

Taconis® has recorded a number of ingenious ideas in 
the patent literature which utilize displacers and regener- 
ators to produce refrigeration at low temperatures. These 
are similar in many respects to the Philips machine and also 
to machines which have been developed recently at 
Arthur D. Little, Inc.7-® 

The most recent developments at Arthur D. Little are 
the results of an effort to overcome the limitations, while 
retaining the most attractive features, of systems utilizing 
thermal regenerators and displacers. 


The single-stage cycle 

The new refrigerator utilizes the feature of reciprocating 
flow through a thermal regenerator rather than steady flow 
through countercurrent heat exchangers. Advantages 
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which can be realized with thermal regenerators as com- 
pared with heat exchangers, such as extreme compactness 
and simplicity of construction, are thereby preserved. 
Unlike the Kirk—Philips cycle, however, the present re- 
frigerator uses a conventional gas compressor with no 
phase co-ordination between the compression and expan- 
sion sections of the system. Another point of distinction is 
that mechanically controlled intake and exhaust valves, 
operating at room temperature, are used to control the 
flow of gas through the refrigerator and to define the pres- 
sure—volume cycle in a manner similar to a conventional 
expansion engine. A unique feature of the new refrigerator 
is that the expansion engine does not deliver mechanical 
work to the exterior; instead, it rejects enthalpy in the form 
of heat. The somewhat remarkable consequence of this is 
that the exhaust stream from the refrigerator is heated to a 
higher temperature than that of the entering stream. 


Intake 
\ valve 


Compressed 
gas supply 

Expanded 

gas return 


Exhaust 
valve 


Warm end 


Cylinder 


Regenerator Displacer 


Cold end 


Refrigeration 
load 


Figure 1. Single-stage refrigerator 


A single-stage refrigerator utilizing the new cycle is 
shown schematically in Figure 1. It consists of a thin wall 
stainless steel cylinder into which is fitted a displacer in the 
form of a fairly closely fitting plastic piston which is free to 
move longitudinally throughout a limited distance. A rod 
is attached to the room temperature end of the displacer 
and is used to control its motion. The other end of the 
displacer—cylinder combination is at the low temperature. 
As the displacer is moved, the warm and cold volumes at 
the two ends of the cylinder are varied, one increasing at 
the expense of the other. However, these two volumes are 
connected together through a thermal regenerator having 
relatively low resistance to gas flow. Consequently, the 
pressure above and below the displacer is always approxi- 
mately equal; and no appreciable force is required to move 
it. 

The thermal regenerator consists of a thin wall stainless 
steel tube filled with circles of fine metal screen which are 
stacked to form a cylindrical structure just equal to the 
internal dimensions of the regenerator shell. The upper end 
of the regenerator is at approximately room temperature 
and the other end is at the low temperature. In spite of the 
large temperature gradient along the length of the regener- 
ator, there is very little heat flow because of the poor ther- 
mal conductivity of the stacked structure. 
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The intake and exhaust valves, which are located at room 
temperature, are operated by cams attached to a crank 
shaft which controls the motion of the displacer. Com- 
pressed gas is supplied to the refrigerator from a more or 
less conventional closed cycle compression system con- 
sisting of a reciprocating two-stage compressor and the 
usual ballast chambers and coolers. 

The cycle may be considered to be made up of four op- 
erations: Pressurization, Intake, Expansion, and Exhaust. 
At the beginning of the pressurization the displacer is at the 
extreme cold end of the cylinder so that the cold volume is 
a minimum and the warm volume is a maximum. The in- 
take valve is opened and compressed gas rushes into the 
warm volume (and also the free volume of the regenerator) 
until the pressure is equal to that of the supply. The sudden 
increase of pressure in the warm volume of course causes 
adiabatic heating and results in a considerable temperature 
rise of the gas. 

With the intake valve remaining open, the displacer is 
moved so as to enlarge the cold volume and diminish the 
warm volume. This displaces the heated gas through the 
regenerator to thecold end, causing it to decrease in volume 
by an amount proportional to the ratio of the absolute tem- 
peratures at the two ends. Since the intake valve is open 
during the operation, the pressure is constant, and ‘more 
gas flows in through the intake valve and mixes with the 
displaced stream entering the warm end of the regenerator. 
At the end of the intake stroke the cold volume has reached 
its maximum value and is charged with cold gas at the full 
intake pressure. The intake valve is closed at the end of the 
intake stroke. 

The expansion operation is carried out by opening the 
exhaust valve slowly over a period of time when the dis- 
placer is essentially motionless. If the expansion occurs too 
abruptly, poor heat transfer will result between the ex- 
hausting gas and the metal matrix of the regenerator. As 
expansion occurs the gas remaining in the cold volume 
falls to a progressively lower temperature. The gas leaving 
the cold volume is partially warmed by heat removal from 
the working temperature (the refrigeration load in Figure 
1) and is subsequently heated all the way to room tempera- 
ture by passage through the regenerator. 

The cycle is completed by moving the displacer so as to 
expel the cold expanded gas remaining in the cold volume 
at the end of the expansion operation. The exhaust valve 
remains open as the displacer is lowered and is closed after 
the exhaust stroke is completed. The initial condition is 
now restored and the whole sequence of operations is 
repeated. 


Thermodynamic considerations 
Each cycle results in an enthalpy change of 


AH = +(p,—p.)V 


for the gas at the warm end of the refrigerator and a cor- 
responding change of AH = —(p,—p,) V at the cold end, 
where p, and p, are the high and low pressures, and V is the 
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displaceable volume. This can be demonstrated by refer- 
ence to Figures 2(a) and 2(b) which show the pressure— 
volume diagram for a complete cycle for the warm end and 
cold end, respectively. The two diagrams are identical 
except for the fact that the p—V area is circumscribed in the 
anticlockwise direction at the warm end, corresponding 
to a reduction of enthalpy. Since all other processes invol- 
ved are constant pressure processes, and since there is no 
mechanical work delivered, it can be concluded that, 
ideally, an amount of heat (p,—p,) V is pumped from the 
low temperature and rejected at room temperature for each 
cycle. 


Pressure 
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Figure 2. Pressure—volume diagrams of (a) warm volume and (b) cold 


volume 


There are many ways in which the machine can deviate 
from ideal behaviour. The principal sources of inefficiency 
can be divided into two groups, thermal inefficiencies and 
pressure-volume inefficiencies. Thermal inefficiency can 
result from inadequate heat transfer or insufficient heat 
storage capacity in the thermal regenerator or at the site of 
the refrigeration load. Heat leakage from the warm to the 
cold ends of the displacer—cylinder assembly or along the 
regenerator also contribute to thermal inefficiency. On the 
other hand, poorly arranged valve timing or restricted flow 
passages may seriously diminish the actual pressure— 
volume area, thereby reducing the enthalpy change per 
cycle. Finally, unnecessarily large gas volumes in the flow 
passages or in the regenerator, valves, or connecting tubes 
result in gas wastage, because the static volume fills with 
gas and discharges each cycle without contributing to 
enthalpy change. 

Generally speaking, minimization of thermal inefficien- 
cies is incompatible with minimization of pressure-volume 
inefficiencies. For example, an ideal thermal regenerator 
would have extremely small flow passages, a large surface 
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area for good heat transfer, and a large volume of metal for 
good heat storage. This is clearly incompatible with the 
requirement of minimum internal gas volume and mini- 
mum restriction to gas flow. In order to achieve maximum 
efficiency for the refrigerator, it is necessary to optimize 
the design in multidimensional space, taking into account 
at least ten independent parameters. These parameters are: 
length and diameter of regenerator, density, specific heat 
and degree of subdivision of metal matrix, length and fre- 
quency of stroke, length and diameter of piston, and maxi- 
mum and minimum pressure. Some of these problems have 
been considered in more detail in previous publications.78 
An actual single-stage refrigerator shows an efficiency 
which is dependent on the temperature span, the efficiency 
decreasing as the span is increased. Several refrigerators 
have been constructed. which produce 65-70 per cent of 
(p,—p:) V in the form of useful refrigeration, as measured 
by an electrical heater, at approximately 100°K, where 
p, and p, are 300 and 75 |b/in? above atmospheric pressure. 


The multi-stage cycle 


It has been noted that the efficiency of a single-stage 
refrigerator decreases as the temperature span increases. 
This suggests the possibility of using several stages ar- 
ranged end-to-end in order to minimize the span covered 
by any single stage and thus permit the attainment of very 
low temperatures with high efficiency. Fortunately, the 
present cycle lends itself very well to this arrangement, as 
shown schematically in Figure 3. The operating cycle in 
this case is exactly the same as for the single-stage machine. 
Refrigeration is produced at three different temperatures: 
80°K, 35°K, and 14°K in the example shown. The 
relative magnitudes of the displaced volumes determine 
the amount of refrigeration produced at each temperature 
and will therefore vary from one application to another. 

The lowest attainable temperature at the cold end of the 
third stage is determined by the performance of the third- 
stage regenerator. The heat capacity per unit volume of all 
metals diminishes very rapidly at temperatures below 20° K 
and their heat storage ability becomes poor. The difficulty 
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Figure 3. Three-stage refrigerator 
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is further aggravated by the fact that the heat capacity per 
unit volume of helium gas increases in inverse proportion 
to the temperature. As a practical matter it has been found 
that 13-14°K may be regarded as a lower limit if reason- 
able efficiency is to be preserved, although it is possible to 
achieve temperatures of 6-8°K if no heat load is applied. 

In the multi-stage cycle, refrigeration produced at inter- 
mediate temperature levels is used for three general 
purposes. Firstly, it is used to make up for the inefficiency 
of the thermal regenerator which feeds gas to that particu- 
lar temperature level. Secondly, it is used to intercept heat 
which is flowing into the colder zones either by conduction 
along the displacers and cylinders or by radiation from the 
warm enclosure. And finally, in the event of the machine 
being used as a liquefier, refrigeration at the intermediate 
temperatures is used to cool down that portion of gas 
which is actually liquefied so as to minimize the refriger- 
ation load on the final stage. 

Heat transfer between the thermal load (regardless of its 
origin) and the refrigeration source is established by means 
of heat transfer stations which form a part of the flow 
channels connecting the regenerators to the expansion 
chambers at the cold end of each cylinder. Such a heat 
exchanger is shown in Figure 1, but they are omitted from 
Figure 3 for simplicity. 


The helium liquefaction cycle 


A three-stage refrigerator such as is shown diagram- 
matically in Figure 3 can be modified by the addition of a 
simple Joule-Thomson cycle to achieve helium lique- 
faction as shown schematically in Figure 4. The Joule— 
Thomson stream, flowing continuously, is cooled by 
countercurrent flow with the unliquefied portion of gas and 
also by bringing it into contact with the heat transfer 
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Figure 4. Helium liquefaction cycle 
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stations of the three-stage refrigerator at 80°K, 35°K, and 
14°K. If helium gas at 300 Ib/in? pressure is cooled to 14°K 
and expanded to | atm in a simple Joule~Thomson cycle, 
approximately 10 per cent of the gas is liquefied; the 
remaining 90 per cent is recirculated and recompressed. 
The overall recirculation system is shown in Figure 5. 
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Figure 5. Helium circulation system 


Table 1. Description of Helium Refrigerator 


Model | Model 2 
Compressor 
Water-cooled Air-cooled 
Joule-Thomson cylinder 
Intake (gauge) 0 Ib/in? 0 Ib/in? 
Exhaust (gauge) 75 |b/in? 75 |b/in? 
Refrigerator cylinder 
Suction (gauge) 75 |b/in? 75 |b/in? 
Exhaust (gauge) 300 Ib/in? 300 Ib/in? 
Volumetric displacement 
J-T cylinder 7-63 ft/min 3-85 ft?/min 


Refrigerator cylinder 15-30 ft?/min 3-85 ft?/min 


Power consumption 10 h.p. 3 h.p. 
Weight: motor and 
compressor 359 Ib 183 lb 
Size: motor and 
compressor 28 in. x 174 in. 20 in. x 124 in. 
x 14 in. xo Dain: 
Refrigerator 
Wattage 
at 14°K 15 W 5 W 
at 4:2°K 4W 0-75 W 
at 3:2°K 3W 0-22 W 
at 2:8°K 0-5 W 
Size of refrigerated space 1 ft® 0-1 ft® 
Size of refrigerator 5 ftx2 ft 27 in. x 15 in. 
x2 ft x 9 in. 
Weight of refrigerator 500 Ib 90 Ib 
Length of 
Cylinder 1 10 in. 7 in. 
Cylinder 2 20 in. 12 in. 
Cylinder 3 30 in. 17 in. 
Diameter of 
Cylinder 1 1? in. 1} in. 
Cylinder 2 1} in. $ in. 
Cylinder 3 fin. +5 in. 
Stroke 1 in. 1 in. 
Revolution rate 75 rev/min 105 rev/min 


See 
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Several closed cycle refrigeration systems have been 
built to date which have operated successfully at temper- 
atures of 4-2°K and lower. While each of these experi- 
mental machines has differed from the others in some 
respects, they have all had certain characteristics in com- 
mon. In each case the system was closed and indifferent to 
atmospheric pressure. The three operating pressures are 
established by pressure regulators which either draw gas 
from or deliver gas to a small storage tank permanently 
connected to the system. In all cases the compressors have 
been standard commercial Freon compressors somewhat 
modified for compression of helium. Special provisions 
have had to be made for adequate lubrication and oil 

separation from the compressed helium. Figures 6, 7, and 


Figure 6. Internal parts of 3 W refrigerator 


8 are photographs of the refrigerator, the refrigerator in 
its enclosure, and the compressor system, respectively, 
referred to as ‘Model 1’ in the preceding Table. Figures 9 
and 10 are photographs of ‘Model 2’ refrigerator and its 
air-cooled compressor system. 

Operating experience to date has indicated that it is 
quite possible to run such refrigerators for periods of 
several hundred hours with little or no attention. Actually, 
the duration of a continuous operation is at present de- 
termined by the efficiency of the oil filters used on the high 
pressure helium supply line. Inadequate filtration of 
entrained oil droplets and oil vapour eventually cause 
clogging of the Joule-Thomson heat exchanger. It is 
probable that another order of magnitude improvement in 
reliability would require the development of compressors 
which do not require liquid lubricants. 


Much of the work reported was sponsored by the 
International Business Machines Corporation and the 
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Figure 7. Complete 3 W refrigerator 


Figure 8. Compressor for 3 W refrigerator 
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Figure 9 (right). 0:75 W refrigerator 


Figure 10 (below). Compressor assembly for 0-75 W refrigerator 


70 
CRYOGENICS + DECEMBER 1960 


Abstracts 


A CRYOSTAT FOR ELECTRICAL 
MEASUREMENTS ON SEMICONDUCTORS 


B. V. ROLLIN J. R. MILLS J. P. RUSSELL 
THE cryostat described is used for routine measurements of the 
resistivity and Hall coefficient of samples of semiconductors. The 
design enables the crysotat to be dismantled and assembled in a 
very short time, so that a different specimen may be measured every 
day. The difficulties of siphoning liquid helium have been avoided by 
using the transport Dewar as the experimental Dewar. 


DENSITY OF LIQUID OXYGEN AS A FUNCTION 
OF PRESSURE AND TEMPERATURE 


A. VAN ITTERBEEK O. VERBEKE 

A DESCRIPTION is given of an apparatus for measuring the density 
of liquefied gases as a function of pressure. Experimental results 
of the density of liquid oxygen as a function of temperature and 
pressure are given from 90° K to 65° K and up to 150 kg/cm?. Some 
thermodynamic properties, such as y, the Poisson’s ratio, and W, 
the velocity of sound, are calculated. A comparison with the existing 
measurements of these thermodynamic functions is made. 


A LOW TEMPERATURE PLANT FOR THE 
PRODUCTION OF HEAVY WATER 


J. HANNY 


THE article describes a plant for the production of heavy water (D,O) 
which has been erected by Sulzer Brothers for Emser Werke AG, at 
Domat/Ems, Switzerland. The raw material is natural water which 
first undergoes seven-fold deuterium enrichment in three electrolysis 
stages. The gaseous mixture of hydrogen and deuterium produced by 
electrolytic dissociation is cooled to — 250°C and thereby liquefied. 
The liquid then goes to a Kuhn rectifying column where the exchange 
processes taking place at 23°K result in 400-fold deuterium enrich- 
ment. The product of the column is burnt with oxygen and is there- 
upon passed to a second Kuhn column in which the heavy water 
concentration is brought up to 99:8 per cent. It is the only known plant 
for the liquefaction of hydrogen which utilizes expansion turbines 
only, without any preliminary nitrogen cooling. 


THE CRITICAL FIELD OF SUPERCONDUCTING 
ALUMINIUM UNDER PRESSURE 


D. GROSS J. L. OLSEN 

THE effect of pressure p on the critical field H, of aluminium has been 
measured between 0-3°K and the transition temperature T¢. 
Pressures up to about 1,700 atm were made by the ice-bomb tech- 
nique. We find 


0H,/ép = —[3:140-2+(1-:7+0-2)(7/T,)*] x 10-® gauss dyn-! cm? 


This yields 
(dIny)/(diInv) = 8+4 


where yT is the electronic specific heat per mole. Results on zinc are 
also reported. Here 0H,/0p appears to be anisotropic. 
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STUDIES OF THE LOW TEMPERATURE 
DISTILLATION OF HYDROGEN ISOTOPES 


T. M. FLYNN 


EXPERIMENTAL data obtained from the operation of 1-35 in. diameter 
distillation columns was used for the design of a pilot plant column. 
This column has been constructed and operated to study the separa- 
tion of hydrogen deuteride from hydrogen by fractional distillation. 

The pilot plant column contains thirty plates 6 in. in diameter 
made of 40 mesh screen, and has a maximum feed rate of 14-2 
ft?(n.t.p.)/min. The column was operated with vapour velocities 
ranging from 0:2 to 2:0 in./sec, and reflux ratios from total to an 
(L/D) of 2:05. Overall plate efficiencies were investigated and found 
to vary, as a function of vapour velocity, from 45 to 55 per cent. 

The overall plate efficiencies of the pilot plant column are sig- 
nificantly lower than those of the smaller diameter columns. This 
difference in plate efficiency is examined theoretically. A model 
considering the effect of liquid mixing on the plate satisfactorily 
describes the change in efficiency and predicts that no further decrease 
in efficiency with increasing diameter will occur. 


COOLING BY ADIABATIC MAGNETIZATION 
OF SUPERCONDUCTORS 


M. YAQUB 


IN order to develop a method of attaining temperatures below 1°K, 
a cryostat has been designed to investigate the fall in temperature on 
magnetizing superconductors. It is found that by magnetizing a 
cylindrical specimen of pure tin from an initial temperature of 0-9°K, 
the lowest temperature which can be reached is 0:35°K. This is sub- 
stantially higher than the expected value of 0:02° K for ideal cooling. 
Magnetizations have also been performed on tin-indium alloys. The 
lowest temperatures obtained with 0:5, 1-0, and 2:0 per cent alloys are 
0:19, 0-18, and 0-345°K, respectively. It has been shown that the time 
taken to warm from 0:35°K, the lowest temperature attained with 
pure tin, can be as long as six hours. The causes of irreversibility, 
supported by indirect evidence, have been discussed. 


SUPERFLUIDITY AND SPECIFIC HEAT OF 
LIQUID HELIUM IN ‘VYCOR’ POROUS GLASS 


D. F. BREWER’ D.C. CHAMPENEY K. MENDELSSOHN 


THE flow rate of liquid helium through disks of ‘ Vycor’ porous glass 
has been measured as a function of temperature at constant pressure 
head. The specific heat of liquid helium within the same disks has 
also been measured. The mean pore diameter as given by adsorption 
measurements was 71 A, and the onset of superfluidity and the 
specific heat maximum occurred, respectively, at 2.0+0-1°K and at 
2:06°K. The effect of pore size distribution is discussed, and the 
results compared with similar measurements using the unsaturated 
helium-II film. 


DISCUSSION ON ENTHALPIES OF MIXTURES 
M. RUHEMANN 


Commission I of the International Institute of Refrigeration, meeting 
at Eindhoven in June, 1960, devoted a session to the lack of infor- 
mation on the enthalpies of gaseous and liquid mixtures, under the 
chairmanship of the author. A short report is given here. 
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Resumes 


CRYOSTAT POUR MESURES ELECTRIQUES 
SUR LES SEMICONDUCTEURS 


B. V. ROLLIN J. R. MILLS J. P. RUSSELL 


Le cryostat décrit est utilisé pour les mesures de la résistivité et du 
coefficient Hall d’échantillons de semiconducteurs. La construction 
permet de démonter et d’assembler le cryostat en un temps tres court, 
ce qui permet de mesurer chaque jour un spécimen différent. Les 
difficultés de transvasement de I"hélium liquide par siphon ont été 
évitées en employant le Dewar de transport comme Dewar expéri- 
mental. 


DENSITE DE L’OXYGENE LIQUIDE EN 
BONCTION, DESLALERESSION ET.ADESLA 
TEMPERATURE 


A. VAN ITTERBEEK O. VERBEKE 


UNE description est donnée d'un appareil permettant de mesurer 
la densité de gaz liquéfiés en fonction de la pression. Des résultats 
d’expérience sont donnés pour la densité de l’oxygene liquide en 
fonction de la température et de la pression entre 90° K et 65° K 
et jusqu’a 150 kg cm?. Certaines propriétés thermodynamiques sont 
calculées, telles que y, le rapport de Poisson, et W, la vitesse du 
son. Une comparaison est faite avec les mesures existantes de ces 
fonctions thermodynamiques. 


UNE INSTALLATION POUR LA PRODUCTION 
D*EAU LOURDE A BASSE TEMPERATURE 


J. HANNY 


L°ARTICLE qui Suit décrit une installation pour la production d’eau 
lourde (D,O). établie par Sulzer Fréres dans les usines d’Ems, a 
Domat Ems (Suisse). La matiére premiére est l'eau naturelle, 
enrichie sept fois en deutérium par trois phases d’électrolyse. Le 
mélange gazeux dhydrogene et de deutérium, provenant de la 
décomposition électrolytique. est refroidi a —250 C et liquéfié, 
puis ce liquide passe dans une colonnea rectifier Kuhn, ou |’opération 
déchange a 23 K enrichit de 400 fois la teneur en deutérium. Le 
produit ainsi obtenu est brilé avec de l‘oxygene puis concentré dans 
une seconde colonne Kuhn, pour laisser en fin de compte de l'eau 
lourde a 99.8 pour cent. C’est a notre connaissance la seule instal- 
lation dans laquelle on n applique que des turbines de détente. sans 
aucun refroidissement préliminaire de |’azote. 


LE CHAMP CRITIQUE DE L’ALUMINIUM 
SUPRACONDUCTEUR SOUS PRESSION 


D. GROSS OLSEN 


L’ErFET d'une pression p sur le champ critique H. de |’aluminium a 
éte mesure entre 0.3 K et la température de transition 7. Des 
pressions allant jusqu’a environ 1.700 atmospheres ont été réalisées 
au moyen de la technique de la ‘bombe a glace’. On constate que 
0H. Op = —[(3,1+0.2+(1.7+0.2)(T/T.)?] x 10 * gauss dyne } cm? 
ce qui donne 

(d In y)(d Inv) = 8+4 


yT etant la chaleur spécifique électronique par mole. On donne 
également les résultats relatifs au zinc. Dans ce cas, 0H../0p semble 
€tre anisotrope. 


ETUDES DE LA DISTILLATION A BASSE 
TEMPERATURE DES ISOTOPES D’HYDROGENE 


l. M. FLYNN 


Lrs données expérimentales obtenues au cours du fonctionnement 
de colonnes de distillation de 1,35 in. de diamétre ont servi a |’étude 
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d'une colonne pour installation pilote. Cette colonne a été con- 
struite et utilisée pour l’étude de la séparation de HD et H, par 
distillation fractionée. 

La colonne de l’installation pilote est constituée de trente plateaux 
de 6 in. de diamétre, faits de toile 4 tamis 4 40 mailles par pouce, 
le taux d’alimentation maximum est de 14,2 ft? par minute. La 
colonne fonctionne avec des vitesses de vapeur allént de 0,2 a 
2,0 in./sec et des taux de reflux variant d°un reflux total a un (L/D) 
de 2,05. L’efficacité moyenne des plateaux a été étudiée et ona 
constaté qu'elle varie, en fonction de la vitesse de la vapeur, de 
45 a 55 pour cent. 

L’efficacité moyenne des plateaux de la colonne de I’installation 
pilote est nettement inférieure a celle des colonnes de diamétres 
plus faibles. L’auteur examine, d'un point de vue théorique, cette 
différence dans l’efficacité des plateaux. Un modéle tenant compte 
de leffet du brassage du fiquide sur les plateaux permet de rendre 
compte de facon satisfaisante du changement d’efficacité et de 
prévoir qu'une augmentation du diamétre ne déterminera pas une 
nouvelle réduction de Iefficacité. 


REFROIDISSEMENT PAR L’AIMANTATION 
ADIABATIQUE DES SUPRACONDUCTEURS 


M. YAQUB 


PouR mettre au point une méthode qui permettra d’atteindre des 
températures inférieures a 1° K, on a concu et construit un cryostat 
pour l’étude de l’abaissement de température qui accompagne 
l'aimantation des supraconducteurs. On a constaté que la tempéra- 
ture la plus basse qu'on puisse atteindre en aimantant un échantillon 
cylindrique d’étain pur a partir d’une température initiale de 0,9° K, 
est de 0,35° K. On espérait atteindre 0,02° K, pour le refroidissement 
idéal, et la température obtenue est évidemment bien plus élevée. On 
a également aimanteé des alliages d’etain et d’indium. Les tempéra- 
tures les plus faibles que l’on ait réussi a atteindre, avec des alliages a 
0,5, 1 et 2 pour cent ont été de 0,19, 0,18 et 0,345° K respectivement. 
On a montré que le temps nécessaire pour le réchauffage, a partir de 
0,35 K, la température la plus basse atteinte avec l’etain pur, peut 
aller jusqu’a six heures. L’auteur examine les causes de l’irré- 
versibilité, en s'appuyant sur des renseignements indirectes. 


SUPERFLUIDITE ET CHALEUR SPECIFIQUE 
DE L’HELIUM LIQUIDE DANS UN VERRE 
POREUX ‘VYCOR’ 


D. F. BREWER D.C. CHAMPENEY K. MENDELSSOHN 


ON a mesuré |’écoulement de I’hélium liquide 4 travers des disques 
de verre poreux *Vycor’ en fonction de la température, sous une 
Surpression constante. On a également mesuré la chaleur spécifique 
de l’hélium liquide dans ces mémes disques. Le diamétre moyen des 
pores, donné par des mesures d’adsorption, était de 71 A; la super- 
fluidité apparaissait a 2,0 + 0,1’ K et la chaleur spécifique passait 
par un maximum a 2,06” K. Les auteurs examinent l’influence de la 
répartition en dimension des pores et comparent les résultats a ceux 
des mesures semblables faites avec un film d’hélium-II non saturé. 


DEBAT SUR LES ENTHALPIES DES MELANGES 
M. RUHEMANN 


LA COMMISSION | de l'Institut International de Réfrigération, réunie 
a Eindhoven, en juin 1960, a consacré toute une séance au manque 
de renseignements sur les enthalpies des mélanges gazeux et 
liquides. L’auteur, qui a présidé cette seance, donne, ici, un bref 
compte rendu de la discussion. 
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Ausztuige 


EIN TIEFTEMPERATURGERAT FUR 
ELEKTRISCHE MESSUNGEN AN HALBLEITERN 


B. V. ROLLIN J. R. MILLS J. P. RUSSELL 


Das hier beschriebene Tieftemperaturgerat dient zu Routinemessun- 
gen des spezifischen Widerstandes und des Hall-Koeffizienten von 
Halbleiterproben. Das Gerat ist so konstruiert, dass es sehr schnell 
auseinandergenommen oder zusammengesetzt werden kann, so dass 
taglich eine andere Probe’gemessen werden kann. Die sich beim 
Uberhebern von Helium ergebenden Schwierigkeiten wurden durch 
Benutzung des Transportdewargefasses als Messgefadss uberwunden. 


DICHTE DES FLUSSIGEN SAUERSTOFFS ALS 
FUNKTION VON DRUCK UND, TEMPERATUR 


A. VAN ITTERBEEK O. VERBEKE 


BESCHREIBUNG eines Gerdtes zur Messung der Dichte verfliissigter 
Gase in Abhangigkeit vom Druck. Die Dichte des fliissigen Sauer- 
stoffs wurde zwischen 90° K und 65° K und bis zu 150 kg/cm? als 
Funktion der Temperatur und des Druckes bestimmt. Einige thermo- 
dynamische Eigenschaften wie z.B. y, das Poissonsche Verhiltnis, 
und W, die Schallgeschwindigkeit, wurden berechnet und mit frii- 
heren Werten dieser thermodynamischen Funktionen verglichen. 


EINE TIEFTEMPERATURANLAGE ZUR 
HERSTELLUNG VON SCHWEREN WASSER 


J. HANNY 


Die folgenden Ausfiihrungen beschreiben eine Anlage zur 
Gewinnung von schweren Wasser, welche in den *‘Emser Werken 
AG’ in Domat/Ems steht und von Gebriider Sulzer hergestellt 
wurde. Als Ausgangsprodukt dient natiirliches Wasser, das sich in 
drei Elektrolyseurstufen siebenfach an Deuterium anreichert. Das bei 
der elektrolytischen Zerlegung entstehende gasformige Wasserstoff- 
Deuterium-Gemisch wird auf — 250°C abgekihlt und verflissigt. 
Anschliessend gelangt die Fliissigkeit in eine Kuhn-Rektifizier- 
kolonne, wo der Austauschvorgang bei 23°K eine vierhundertfache 
Deuteriumanreicherung bewirkt. Das gewonnene Produkt wird mit 
Sauerstoff verbrannt und in einer zweiten Kuhn-Kolonne zu 99,8 
prozentigem schwerem Wasser aufkonzentriert. Es ist die einzige 
heute bekannte Anlage zur Verfliissigung von Wasserstoff, welche 
nur Expansionsturbinen und keine Stickstoffvorkihlung bendtigt. 


DIE KRITISCHE FELDSTARKE VON UNTER 
DRUCK STEHENDEN ALUMINIUM- 
SUPRALEITERN 


D. GROSS J. L. OLSEN 


Der Effekt des Drucks p auf die kritische Feldstarke H. von Alumin- 
ium wurde zwischen 0,3° K und der Ubergangstemperatur. 7. 
gemessen. Driicke bis zu 1700 atm. wurden mittels der Eisbomben- 
methode hergestellt. Es ergab sich 


0H./op = —[3,1+0,2+(1,7+0,2) (7/T.)?] x 10 * gauss dyn ! cm? 


und daraus 
(d In y)/(d In v) = 8+4 


worin yT die spezifische Elektronenwarme pro Mol ist. Ergebnisse 
mit Zink werden ebenfalls besprochen, wobei @H,/dp anisotropisch 
zu sein scheint. 
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TIEFTEMPERATURDESTILLATION VON 
WASSERSTOFFISOTOPEN 
T. M. FLYNN 


VERSUCHSWERTE, die im Betrieb mit einer Destillationssaule von 
1,35 in. Durchmesser erzielt wurden, sind fiir die Bemessung einer 
Versuchkolonne verwertet worden. Diese Kolonne wurde gebaut 
und betrieben, um die Trennung von Wasserstoffdeuterid und 
Wasserstoff durch fraktionierte Destillation zu studieren. 

Die Versuchsanlage enthdlt dreissig Boden von 6 in. Durch- 
messer aus 40-er Maschensieb und hat eine maximale Zuflussmenge 
von 14,2 ft?/min. Die Kolonne wurde mit Dampfgeschwindigkeiten 
von 0,2 bis 2,0 in./sec und mit Ruckflussverhdltnissen von unendlich 
bis zu (F/E) von 2,05 betriecben. Messungen des Verstarkungs- 
verhaltnisses ergaben Schwankungen als Funktion der Dampf- 
geschwindigkeit zwischen 45 und 55°,. 

Das Verstaérkungsverhaltnis der Versuchsanlage was bedeutend 
niedriger als dasjenige von Kolonnen mit kleinerem Durchmesser. 
Dieser Unterschied im Verstaérkungsverhaltnis wird theoretisch 
untersucht. Ein die Wirkung der Flissigkeitsmischung auf dem 
Boden veranschaulichendes Modell beschreibt die zufriedenstellend 
die Anderung des Verstarkungsverhdltnisses und ldsst darauf 
schliessen, dass bei noch grdsseren Durchmessern keine weitere 
Abnahme stattfinden wird. 


KUHLUNG DURCH ADIABATISCHE 
MAGNETISIERUNG VON SUPRALEITERN 
M. YAQUB 


Um ein Verfahrens zur Erreichung von Temparaturen unter | K 
zu entwickeln wurde ein K4lteregler entworfen, der den Tempera- 
turabfall bei der Magnetisierung von Supraleitern zu studieren 
gestattet. Es wurde gefunden, dass durch Magnetisierung eines 
Zylinders aus reinem Zinn bei einer Ausgangstemperatur von 0,9 K 
eine tiefste Temperatur von 0,35 K erreicht werden kann. Diese ist 
bedeutend hoher, als der bei idealer Kiihlung erwartete Wert von 
0,02 K. Magnetisierungsversuche wurden auch mit Zinn-Indium- 
Legierungen .durchgefiihrt. Die tiefsten Temperaturen, die mit 
0,5°,, 1,0°5 und 2,0°,, Legierungen erzielt wurden, sind 0,19, 0,18, 
bzw. 0,345’ K. Es hat sich gezeigt, dass die Aufwarmzeit von 
0,35° K, d.h., von der tiefsten mit reinem Zinn erreichten Tempera- 
tur, bis zu sechs Stunden betragen kann. Die Ursachen der [rreversi- 
bilitat, die durch indirekte Beweise gestitzt sind, wurden besprochen. 


SUPRAFLUIDITAT UND SPEZIFISCHE WARME 
VON FLUSSIGEM HELIUM IN POROSEM 
‘“VYCOR’ GLAS 
D. F. BREWER D. C. CHAMPENEY K. MENDELSSOHN 


Die Durchflussleistung von fliissigem Helium durch pordése * Vycor’ 
Glas-Platten wurde als Funktion der Temperatur bei konstanter 
Druckhohe gemessen. Die spezifische Warme von fliissigen Helium 
innerhalb derselben Platten wurde ebenfalls gemessen. Der durch 
Adsorptionsmessungen bestimmte durchschnittliche Porendurch- 
messer ergab sich zu 71 AE. Der Beginn der Suprafluiditat, bzw. 
das Maximum der spezifischen Warme, traten bei 2,0 + 0,1 K und 
bei 2,06° K ein. Der Einfluss der Verteilung der Porengrossen wird 
besprochen, und das Resultat wird mit 4hnlichen Messungen ver- 
glichen, die mit ungesattigtem Helium-II-Film erhalten wurden. 


AUSSPRACHE UBER DIE ENTHALPIEN VON 
GEMISCHEN 
M. RUHEMANN 
Die |. Kommission des Internationalen Kalte-Instituts hat sich in 
ihrer Sitzung in Eindhoven im Juni 1960 unter dem Vorsitz des 
Verfassers mit den unzureichenden Kenntnissen der Enthalpien 
von Gas- und Flissigkeitsgemischen befasst. Der Verfasser erstattet 
einen kurzen Bericht. 
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Pe3iome 


KPHOCTAT JIA dJIEKTPMUECKUX A3MEPEHMA 
MOJYNWPOBOHAKOB 


B. V. ROLLIN J. R. MILLS J. P. RUSSELL 


OnvcaHHHit KpuoctaT ynoTpeOmAeTcCAH JIA O4epeAHbIX M3- 
MepeHMit yAeubHOrO COMpoTMBIeEHMA M KOappuuMenta NOMIAa 
o6pa3noB nomympoBpoyHHKoB. baarojapaA ero KOHCTpyKUMK 
KpMocTaT MO?KHO pas00OpaTb uM cOOpaTb B TeYeHHe OYeHb 
KOPOTKOrO BpeMeHMM, TAK UTO e*KeHEBHO MO?KHO H3MePATb 
pasHpe o6pasubt. ITlomb30BaHMeM TpaHCNOPTHHIM {l0apoM 
BMeCTO JKCHEPHMMeHTAIbHOrO Wioapa MO*KHO H30e#HATb TPY.1- 
HOCTH CHPOHMpOBAaHMA PKUAKOTO reaA. 


IJIOTHOCTb ¢HUJKOTO KUC.IOPO.LA RAK 
®VHHKUUA JABJIEHWVA WU TEMMEPATY Pbl 


A. VAN ITTERBEEK O. VERBEKE 


jlaetcA omuMcaHMe annapaTa (JIA W3MepeHHA MAOTHOCTM 
Pas¢KWVKEHHEIX Ta30B Kak QMYHKIMA WaBeHHA. lawTca 
JKCHEPMMeHTAJIbHBIe Pe3yIbTAThI O MIOTHOCTH *KILIKOFO KMCIO- 
pogwa Kak P@yHKUMA TemMepaTyph Mm ZaBnenua oT 90° K yo 
65° K u go 150 xr/cm?. Busruncaenbl HeKOTOpble TepMOsM- 
HaMMYeCKMe CBOMCTBa, KaK HanpuMep y, oTHOWeHHeE Ilyacona 
u W, ckopocTb 3ByKa. OTM TepMOAMHaMMYecKMe PyHKUMM 
CPpaBHUBawTCA C CYUIECTBYWWIMMM U3MepeHuAMM. 


HU3SKOTEMITEPATYPHAA YCTAHOBKHA JIJIA 
MPOU3BO]CTBA TAKEJION BOLI 


J. HANNY 


B cratbe OMMChIBaeTCA YCTAHOBKA AIA NpPOUSBOACTBA TAKEO 
Bogbt (D2gO), KoTopan ycraHoBneHa d@upmow Cyauep (Sulzer 
Brothers) qua @upmbr Omcep Bepxe AI’ (Emser Werke AG) 
y Jjomar/Omc B Ipeittuapuu. Crrpbem cayKuT npMpoqHaA 
Boa, KOTOpaA cnepBa nosBepraeTcH CeMMKpaTHOMy OOora- 
I€HMIO OTHOCMTeEAbHO eliTepuA B TpéX CTYMeHAX 9TeKTpO- 
mm3a. Ta3s0BpaH cmMecbh BOWOpOAa-elitepuA, mpousBesqeHHaA 
JAEKTPOIMTUYECKON TUCcouMalMel, oxanKqaerca Wo — 250° C 
u Tak omMHaeTCA. IIOTOM *UKOCTh IPUBOAMTCA B pPeKTUU- 
KalMOHHYIO KONWOHY Tuna Hy, B KOTOpOw ocTaBaeTcA 
400-KpatHoe oO6oraujeHue OTHOCMTeIbHOrO felTepuA myTém 
OOMeHHEIX Mpoyeccos npu 23° K. I[lpoxynr o6pasopanuniit 
B KOWIORe CKMTAETCA IPM KUCTOPOAe MU 3aTeM NPMBOAUTCA B 
apyryio KOomHOHY Tuna Kyu, B KOTOpOw KOHWeHTpalMA 
TAKeION BOLI yBemmunuBaeTtcA JO 99,8%. Onucripaeman 
CTAHUMA ABIACTCA CE MHCTBEHHOM YCTAHOBKOM WIA OmKMrKeHVA 
BOJOpOLa Ha KOTOPOM MpMMeHAIOTCA TOIbKO TYPOMHEI C pac- 
wMpeHveM napa 6e3 npesqBapMTeMbHOTO OXTAKTeCHMA asOTOM, 


KPUTMYUECHKOE TOJIE CBEPXMPOBOJJAIN ETO 
AJIIOMMHMA TO TABJIEHUEM 


D. GROSS J. L., OLSEN 


BuIMAHMe JaBIeHMA p Ha KpuTMuecKoe none He anloMMHNA 
OBO OpeereHO NpM TemMepatypax mMexKy 0,3° KH um nepe- 
xOMHOH Temnepatypowt T.. Jleqanaa 6Oom6a Onia mpn- 
MeHeHA JIA CO3qaHMA JaBnenni 701.700 ar. Haxogum 


0H,/¢p = —[38,1+0,24+(1,7+0,2) (T/T¢)2] x 10-9 
rayce fwHa-! cm2 
OTO FaerT 
(dIny)/(dInv) = 8+4 
rye y7'— slekTpOHHaA YilelbHaH Tellofa Ha rpaMM- 
MONeKYIY. PesynbTaThl IA WAHKA TOME POKIAAbIBAOTCA. 
B oToM 2ke Cilyyae 0H e/Ep KamkeTCA aHM30TPONHLIM. 
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U3Y4YEHMA HU3SKOTEMMEPATYPHONM JMCTUUI- 
JIAUMUMA BOJOPOJHbIX M380TOIOB 
T. H. FLYNN 

OnsITHble aHHble, MOMy4YeHHble OT paOOTh AMCTMIIAWMOHHBIX 
KONOH MamMeTpom 1,35 joliMa, MpMMeHMIMCbh JIA KOH- 
CTPYKUMM NOys8aBOACKOM KONIOHBI. ITA KOIIOHA MOCTpoeHa 
If IpMBeyena B MeMCTBMe AIA U3y4eHMA BbIAeLeHAA neuTepuya 
BO_OpoAa U3 BOMOpoAa nyTeM PpakUMOHMpOBaHHOH AMCTUI- 
AAW. 

TlonysapogckKaA KONOHA COMepKUT TpUAWaTb TapeOK 
auametpom 6 juoiimMoB, CflewaHHEIx “3 cur 40 orBepcernit; 
WakCHMasIbHa cKOpocts nuTanua — 14,2 dyron® (H.T.2.) B 
winHyTy. Konnona padorama mpa cKopocrax mapa oT 0,2 70 
2,0 qwoiima B CeKYHAY mM mpm uncnax lerMbl (+ KOCTH / 
anetnaaatT) oT nomHoTH go 2,05. OOdume aphexruBxocTu 
TapetOk ObLIM MCCAeOBaHbl, HAXOAMIOCh UTO OHM U3MCHAWT- 
cA, Kak PyYHKUMA CKOpocTH napa, oT 45 Wo 55%. 

O6ume KoadiuywenThl mouesHoro jelCTBMA TapeJOK 
NOWy3aBOACKOH KOMNOHBI 3HAYMTebHO HMKe 4eM K.II.]L. 
TapesOK KOJINOH MeHbUIMX [MaMeTpOB. OTA pasHMija Mey 
K.IL.J,. Tapesiok MecrefqyetcA Teopetuyecku. Mogemb, 0 
KOTOPOH YUUTHBAeTCA BIMAHMe CMeIIMBAaHMA *KUKOCTH Ha 
TapesIKy, MOKa3bIBaeT M3MeHeHHE HK... yOBMeTBOPUTeILHO 
M MpewcKasblBaeT 4TO OOAbWerO YMeHbIEHHMA K.1.4. He OyAeT 
pv yBesIm4eHuu Tuametpa. 


OXJIAMJEHWUE 10 METOJY AJMABATUYUECKON 

MArHETHM3AUUU CBEPXIPOBOTHYUKOB 
M. YAQUB ‘ 
UToObl pasBuBaTb MeTOX MOCTWKeHMA TemMMepaTyp HuUrKe 
1° HK, cKOHCTpyMpoBaH KpMOCTaT JIA WCCMeqOBaHUA CHMKe- 
HMA TeMMepaTypbl Mp MarHeTU3alMM CBepXMpOBOAHUKOB. 
Haxogutca 470 Ip MarHeTM3alMu WMIMHApMyeckoro OOpasia 
YMcTOrO O0Ba OT UCXOqHOM TemnepatypH 0,9° H Hucuiaa 
qocnHraemaaA TemMueparypa — 0,35° KH; ora remnmeparypa 
CYUIeCTBEHHO BBIMe OKMAaHHOH TemnepatypH 0,02° KH npn 
ueasbHOM OxJaHKeHunM. Buino0HeHbI TOKE MarHeTU3alun 
cnmaBoB On0B0-MHAMit. IIpu cnmaBax 0,5%, 1,0% u 2,0% 
COOTBETCTBEHHbIe HUCHIME JOCTHPHYThIe TeMiepatypb — 0,19° 
H, 0,18° Hu 0,345° HK. Tlokaspipaetca uTo cpok HarpeBaHMA 
ot 0,35° K, Hucmieli TemnepaTypbl AOCTMrHyTOM NpM 4uMcTorO 
O0Ba, MOET YAIMHMTECA JO WecTH yacoB. Obdcymaerca 
NIpM4MHbE HeEOOPpaTMMOCTH, C MOMONIbIO HeMmpAMOTO jOKasa- 
TeIbCTBA. 


CBEPXTEHYYECTS U TENJIOEMKOCTS HU JTKOLO 
TEJIMA B WOPUCTOM CTEHJIE «BUHKOP» 

D. F. BREWER D.C. CHAMPENEY K. MENDELSSOHN 
VMsmepiiacbh CKOpOCTh TeYeHHA *KUAKOLO TeIMA Yepe3 TMCKI 
nmopucToro crekia «BuKOp» Kak PYHKUMA TemMepaTyph mpi 
NOCTOAHHOM Hanope. Viamepumacb TaKeKe TerJOCMKOCTb 
AMUAKOLO PeIMA B TeX Ke CAMBIX AMCKAaX. CpesHuit Juamerp 
nop, Onpexeuénublii M3MepeHHAMM afcopOuun, cocraBnAd 
71 A, uw Hayano cBepxTeKy¥ecTM M MaKCMMaJIbHaA Tens0eM- 
KOCTh HaxogqMIMCch upw 2,0+0,1° K uw npu 2,06° K coorser- 
cTBeHHO. OOcyaeTCA BIMAHMe paciipeqeweHuA pasmepoB 
Nop, M pesyIbTaTbl CpaBHMBAWTCA C MOAOOHKIMM M3MepeHUA- 


MM, MOUYYCHHBIMM NYTEM TPMMeHEHMA HEHACHILCHHOM MIeHKIL 
reamna—lII. 


OBCY#K]TEHWE JHTAJIBIMA CMECEN 
M. RUHEMANN 
Romuccua 1 Mempyxapoguoro Wucruryta no OxsaKennn, 
COCTOABINAACA B OiiHAroBene (B Tommanjun) B uioHe 1960 r., 
oTAada CeccMio OOCYKAeEHMIO BON|POCa HEWOCTaTKa MHPOpMaLMit 
N10 NOBOAY dHTaIbuMit PasOBblX M KMAKMX CMecel, NOW mpey- 
CeqaTeIbCTBOM aBTOpa. Bech NpeACTAaBIEH KpaTKMit OKA. 
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A Cryostat for Electrical Measurements on 
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A FUNDAMENTAL factor in research on semiconductors is a 
knowledge of the concentrations of acceptor and donor 
impurities, and of the activation energy associated with 
these. The usual way of discovering this is to investigate 
the variation of the resistivity and the Hall coefficient 
with temperature between liquid helium temperature and 
room temperature. The cryostat described enables 
temperatures between 4-2° K and 1-5° K to be obtained 
by lowering the vapour pressure above the helium bath, 
and temperatures above 4-2° K to be obtained by evacuat- 
ing the space round the vessel containing the specimen, 
and passing current through a heater. 

One of the major difficulties in the use of liquid helium 
is the transfer of the liquid from one Dewar to another, 
particularly since the cost of helium makes it necessary 
to save all the evaporated gas. In most installations two 
transfer operations are necessary; first from the storage 
container into a transport Dewar and then from the 
transport Dewar into the cryostat. In the present design 
the second operation has been eliminated by enabling 
the transport Dewar to be put on to the cryostat and to 
be used as the experimental Dewar. Loss of gas during 
this operation is prevented by a sliding seal (clearance, 
0-001 in.) between a brass sleeve A and a long thin- 
walled german silver tube B forming the outer casing. 

As materials such as germanium and silicon may have 
a resistivity greater than 10!° Qcm at low temperatures, 
special precautions have been taken in insulating the 
leads, in the electrical screening, and in avoidance of 
frictional charges which can be produced by the boiling 
liquid helium shaking the polythene sleeving encasing 
the leads. This is prevented by enclosing the leads in the 
tube B so that they are not in contact with the liquid 
helium. The specimens may be sensitive to long wave- 
length infra-red radiation from the warmer parts of the 
apparatus and must be shielded. Surface contamination 
and physical strains must be reduced to a minimum. 

The specimen is a bar 1-5 cmx0-3 cmx0-1 cm cut 
from a single crystal ingot, with six contacts enabling 
the resistivity and Hall coefficient to be measured using a 
transverse magnetic field. This specimen is suspended 
from the mount by its electrical leads only, so that there 
is no straining or contamination. The mount is based on 
a modified B7G valve-holder which fits a standard 8-pin 
metal—glass seal, opaque to infra-red radiation, forming 
the top of a small thin-walled brass can C containing the 
specimen. The can C, which can be dismantled ata Wood’s 
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metal joint, is filled with helium gas and is connected toa 
manometer by | mm german silver capillary, thus acting 
as a simple gas thermometer. A carbon resistance 
thermometer is fixed to the bottom of the can C and is 
used as a secondary thermometer. 

The long german silver tube B forms a vacuum space 
in which the can C is suspended by thin german silver 
supports. Heat leakage down the leads has been reduced 
by cementing them to Perspex spacers which are attached 
to copper rings making thermal contact with B by means 
of the springs D. Copper shields reduce the thermal 
radiation down the tube. A piece of 35 s.w.g. constantan 
wire attached to the bottom of the can C presses against 
the bottom of B forming a small heat leak enabling the 
temperature of C to be controlled more readily. At the 


onal 


A 


Figure I. The cryostat 
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top of the cryostat are two standard 7-pin metal-glass 
seals for the electrical connections. 

The liquid helium is contained in a silvered glass 
Dewar E of normal design, surrounded by liquid oxygen 
in a second Dewar F, both having narrow tails to fit an 
electromagnet. The inner Dewar is closed by a brass 
top G (made gas-tight by rubber), to which the sleeve A 
is attached. Evaporated helium gas escapes through the 
tube H and is collected in a balloon pending return to the 
liquefier. Radiation downwards into the liquid helium is 
reduced by the silver-plated brass shield J which is cooled 
by the evaporated gas. A spring skirt K round the edge 
of the shield forces most of the gas to flow between the 
tubes L and B, assisting in the cooling of both B and J 
which reach a steady temperature of about 120° K. 


Procedure 


The cryostat is first of all assembled: the specimen is 
soldered on to the mount which is plugged on to the 
8-pin seal, the can C is soldered together, the tube B is 
pushed up over the springs D and is sealed by a Wood’s 
metal joint. The air is now pumped out of both the can C 
and the tube B; C is filled with helium to a pressure of 
about 30 cm, B with about | mm helium exchange gas. 

The inner Dewar is pre-cooled with liquid oxygen 
(the liquid oxygen being carefully emptied out again), 
and put into the outer Dewar which is then filled with 
liquid oxygen. The required amount of liquid helium is 
then siphoned from the storage vessel into the inner 
Dewar. The Dewars are now positioned under the cryo- 
Stat, a flexible return line is attached, and the Dewars are 
raised until the bottom of the tube B is just inside the 
tube L. The evaporating gas cools the specimen to about 
150° K in 5 min. The Dewars are now lifted into their 
final position. It is very important that this should be 
done slowly as the tube B is then cooled efficiently by the 
evaporating gas and a minimum of liquid helium is used 
in cooling the apparatus to its final temperature. With 1 1. 
of liquid helium the Dewars should be raised 1 in. every 
minute, and about 70 cm® of liquid helium evaporate 
by the time the Dewars are in the final position. When in 
the final position a support on the sleeve A is put into a 
clamp, the outer Dewar is attached to a supporting ring 
by springs, a piece of rubber is slipped over the top of the 
sleeve A and wired on to make it vacuum-tight, and the 
flexible return line is replaced by a metal pipe suitable 
for pumping. 


Associated equipment 


The cryostat has the usual high vacuum system, and a 
pump for lowering the vapour pressure of the helium 
bath. The electrical measurements are made with a 
vibrating reed electrometer in conjunction with a simple 
switching system enabling the voltage across the two 
potential contacts, the Hall contacts, or a comparable 
series resistance to be observed. The system is insulated 
throughout with polythene or polystyrene. Specimen 
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resistances up to 4x 1044 Q can be measured and at this 
resistance Hall voltages of the order of millivolts can be 
observed. 

The carbon resistance thermometer with a heat dis- 
sipation of less than 10-7 W forms one arm of a Wheat- 
stone bridge provided with compensating leads running 
to the bottom of the cryostat. A galvanometer amplifier 
is used to indicate balance and the d.c. output of this 
amplifier is further amplified by a d.c. transistor amplifier. 
The output of this amplifier is passed through the heater 
and automatically controls the temperature. This 
arrangement enables the temperature to be held steady 
to within + 0-2 per cent, which is very necessary because 
of the rapid variation with temperature of the properties 
being observed. 


Performance of apparatus 


Time required for assembly and initial 


cooling 1-1-5 hr 
Volume of liquid helium evaporated dur- 

ing initial cooling 70 cm? 
Evaporation rate of liquid helium 1 cm?/min 


Average duration of experiment using | |. 
of liquid helium 10 hr 
Accuracy of temperature measurement: 
Gas thermometer 
Resistance thermometer 
Temperature stability using amplifier de- 
scribed +0-2 per cent 
Minimum temperature obtainable BSE 


per cent 


+2 
+5. per cent 


(Using liquid hydrogen) (10-8° K) 
Maximum temperature obtainable 290° K 
Leakage resistance 103 © 
Maximum resistance observable 4x 10U QO 
Minimum Hall voltage observable 

At low resistance 0-2 mV 

At high resistance 2 mV 
Principal dimensions 
Length of tube B 80 cm 


Diameter of tube B 3cm 
Wall thickness of tube B 0:4 mm 
Maximum capacity of inner Dewar 1:75 1. 
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IN order to obtain specific heats from the data of sound 
velocity in liquefied gases under high pressure, we have to 
make use of the compressibility coefficient of the liquid. 
However, except for helium and hydrogen, data of this 
kind are rather rare. Because of this, we have started an 
investigation on the density variation of liquefied gases 
such as oxygen, nitrogen, argon, neon, etc., as a function 
of temperature up to 150 kg/cm?. So far measurements on 
liquid oxygen between 90° K and 65° K have been carried 
out. To a first approximation, a straight line is found for 
p=/f(P). From these measurements we have calculated 
the variation as a function of temperature of the isothermal 
(0p/OP)-. 


Experimental method and apparatus 


The apparatus is roughly the same as that which we used 
for our preliminary measurements. It is shown in Figure 1. 


Figure 1. Apparatus for the measurement of density as a function of 


pressure 


We used several calibrated volumes: 


V, the condensation volume. At 273°K and | kg/cm?, V is 
97-17 cm. 

V.., the metal expansion volume. About 41. 

A,, A,, B, C, four glass expansion volumes of different 
sizes. 


Each of these volumes is thermostatted to within 0:01°K. 
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Figure 2. Electrical circuit of the apparatus 


The volume V has been calibrated separately with 
mercury and with distilled water. The dependence of this 
volume on temperature and pressure is calculated with 
the help of experimental data! and of the theory of 
elasticity.2 The pure gas is condensed in this carefully 
calibrated volume V at the boiling temperature. In order 
to determine the mass of the liquid contained in this 
volume, the vessel is warmed up to room temperature and 
the gas expands into V,. With the aid of the P,p,T rela- 
tions at room temperature, it is possible to calculate this 
mass. 

Several corrections have to be made, e.g. gas and liquid 
correction in the capillary tubes. The liquid present in the 
capillary tubes leading to the condensation vessel adds to 
the liquid volume, and this correction is important even 
at low pressure. The temperature of the top of this liquid 
level in the capillary tube is determined from the pressure. 
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Thus an exact knowledge of the temperature gradient 
along the tube is necessary to estimate the height of this 
liquid level. For this purpose we used a set of four thermo- 
couples placed along the tube (Figure 2). The density is 
then given by 
M—-—M,—M, 
poe ORE (ET) 


where M, and M, are the gas and liquid corrections, 
respectively. 

To measure the variation of the density as a function of 
pressure and especially (0p/0P),, we have preferred another 
procedure. Starting with the liquid under high pressure, 
the pressure is successively lowered by expanding a 
quantity of gas into the glass expansion volumes. A de- 
termination of this amount allows calculation of the loss 
of mass and the density variation for each step. When 
vapour tension is reached, a direct measurement with V 
gives the absolute value of the density. 

The temperature is measured with a platinum resistance 
thermometer calibrated with a standard thermometer 
provided by Dr. H. Van Dijk of Leiden, and calibrated 
by him. High pressure measurements are carried out with 
two Heise Bourdon manometers. The oxygen (purity of 
the gas is between 98-9 and 99-2 per cent) used is delivered 
in high pressure cylinders by the firm of Sogaz at Haren 
near Brussels. This gas is condensed twice, so the purity in 
the liquid phase is near 99-95 per cent. 


Table 1. Experimental Results 


4p P. 


p 
(°K) (kg/cm?) (g/cm’) 
64:66 148-17 1:2778 
(Resistance thermometer 136-77 1:2758 
0:760062 2) 125-90 1:2743 
115-32 1:2728 
104-76 1:2701 
83-35 1:2685 
72:57 1:2671 
61-80 1:2657 
50:90 1-2642 
40-12 1:2633 
29-04 1:2621 
7:96 1:2608 
7:82 1:2591 
74:77 “148-92 1:2345 
(Resistance thermometer 138-00 1-2333 
0:99062 Q) 128-00 1-2316 
117-85 1:2300 
107:57 1-2284 
97:29 1:2267 
86:90 1:2249 
76:27 1-2232 
66:72 1-2209 
55:27 1:2197 
44:75 1:2180 
34:20 1-2168 
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Table !—continued 


T - p 
(°K) (kg/cm*) (g/cm*) 
80-18 148-17 1:2138 
(Resistance thermometer 138-02 1:2120 
1:11546 Q) 128-23 1-:2100 
118-12 1:2083 
107-85 1-2063 
97:77 1:2046 
87:53 1:2028 
77:49 1:2011 
67:62 1:1993 
57°53 1:1976 
47-42 11958 
16-35 1-1899 
6:60 1-1883 
85-41 149-10 1-1923 
(Resistance thermometer 139-07 1-:1900 
1:23725 Q) 128-34 1:1879 
117:80 1:1858 
107-10 1:1836 
96:65 1-1817 
86-07 1-1794 
75:46 1:1775 
64-75 1-:1755 
54:20 1-:1734 
43-88 1:1717 
33-20 1:1697 
22:44 1-1680 
12-50 1:1662 
90-29 146-51 1:1707 
(Resistance thermometer 139-95 1-:1665 
1:35252 Q) 125-78 1:1664 
115-87 1-1643 
106-08 1-:1628 
96:10 1-1602 
86-45 1:1582 
76°55 1-1561 
66:93 1-1540 
57:21 1-1520 
47-48 11498 
37:70 1:1483 
27:84 1:1462 
18-18 1:1444 
8-57 1:1424 
1:9] 1:1413 


Discussion and calculation 


The pressure domain is too small for an exact determin- 
ation of a variation of (0p/AP), with pressure. However, 
from our experimental data it is possible to calculate an 
equation of the following form: 


p = A+BP 
with Pin units of kg/cm? and p in g/cm’. Thus we obtained: 
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T (°K) i 

90-29 11-1408 + 1-991 x 10-*x P 
85-41 11634 + 1-904 x 10-4 x P 
80:18 — 1-18724+ 1-787 x 10** x P 
74-77 11-2111 41-589 x 10-4 x P 
64-66  1-25814+ 1-281 x 10-4x P 


We have also calculated, starting from these equations, 
the variation of A and B as a function of temperature 


ai2O 
Ee 
S 
) 
1:27 
A 
1:25 
Q 
1-23 
1-21 
0 50 100 150 
IP ——— (kg /cm?) 


122) 


0 50 100 ISO 


ie SS (kg /em’) 
Figure 3. Experimental results: density of liquid oxygen as a function 
of pressure 
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1-5656 — 4-896 x 10-8 T+.2:136 x 10-* T? 
(—0-53327 +0-02824 T) 10-4 


Il 


From our experimental data the following thermo- 
dynamic properties were calculated: 


W: the velocity of sound; 
Cy, Cp: the specific heats; 
y: Poisson’s ratio. 


This is done with the help of three well-known equations: 


aes ae 2 T (ap 
| SS een (ell 
(ae): p oy] +alari) ee 


This gives Cp(P) and we used known values of the vapour 
pressure curve? as integration constants. Further, we have 


ie eed 3 9) 
"= @plaP); vale) 
and, finally, 
2 
[ssr) ] *? 
aca Be Lara A 3) 
(3p) .*C* 


The results which were obtained are given in Table 2. 
Tables 3, 4, and 5 give the comparisons of our measure- 
ments and other known results. 


Table 2 
T | P| @p/aP nx 10" | (a es mare 
(°K)| (kg/cem?)| (g/cm dyn) (see ref. 3) | (7/See) | yA 
65 0 1-328 0-397 Lig | 1-66 
50 0-396 | 
100 0-394 
150 0-393 
10 0 1472 | 0399 1,060 | 1-66 
50 0-397 | 
100 0-394 
150 0-393 
80 0 1-760 0402 | o71 | 1-66 
50 0-400 | 
100 0-398 
150 0-396 
90 0 2-048 0-406 | 904 | 1-67 
50 0-403 | 
100 0-400. : 
150 0398 ! 
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Table 3. Comparison of Results for Density at the Vapour Tension 


T (? (e 

° (g/cm?) (g/cm) References 
OE) Louvain measurements | Existing measurements 
92-00 1-1335 1-1240 4 
89-00 1-:1470 1-:1479 5 
90-65 1-1395 | 1:1181 6 
90-18 11417 1-1447 i 
91-00 1-1380 1-1415 8 
89-55 1-1448 1:1320 9 
78:15 1-1961 1-:1953 10 


Table 4. Comparison of Results for the Velocity of Sound 


ha Ww 
T P (m/sec) 
(°K) (mm Hg) exp. ve 


(reference 11) 


90:30 771 902:91 902 


77-40 157 1,007-27 See 


- * This large deviation is caused by the fact that Cp is not well known. 
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Table 5. Comparison of Results for Expansion Coefficient 


(aie | (2) 
z v\aT = V\0T Reference 
Our measurements Other authors 


ONNMKNF 


is 


90-31 3:95 4-00 x 10-3 12 
4:12x 107% 5 

77:35 3°81 B85 001 Ons 5) 

90-00 3°95 = 

80-00 3-84 = 

70-00 3-73 = 

65:00 3-68 = 
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HEAVY water is the oxide of deuterium, a hydrogen isotope 
whose nucleus contains a neutron in addition to the normal 
proton, so that its atomic weight is not 1, as with hydrogen 
proper, but 2. Water in which normal hydrogen atoms are 
replaced by deuterium atoms is about 10 per cent heavier 
than ordinary water. In appearance and in chemical be- 
haviour this heavy water can hardly be distinguished from 
H,O, but it nevertheless has certain special virtues and in 
particular is an excellent moderator for nuclear reactors. 

The purpose of a moderator is to slow down by elastic 
impacts the speed of the fast neutrons produced by the 
fission of a nucleus of uranium-235. Since slower neutrons 
are more easily captured by other uranium nuclei, the 
moderator contributes to the initiation and maintenance 
of a chain reaction. It is important, however, that the 
moderator should only slow down the neutrons, and 
should actually capture as few as possible. Heavy water 
captures only one neutron for every five hundred absorbed 
by ordinary water, and the critical quantity of natural 
uranium—the quantity needed to initiate a chain reaction 
—is only 5 tons with a heavy water moderator as compared 
with 50 tons with a moderator consisting entirely of 
graphite. 

If the fuel elements used are enriched with the fissile 
uranium-235, natural water can be employed as a moder- 
ator. Even in reactors of this type, however, the use of 
heavy water is still preferable, since it economizes the 
neutrons on which the chain reaction depends, and 
neutron economy is equivalent to better utilization of the 
available fuel. . 

The principal properties of heavy and natural water are 
compared in Table 1. 


Table 1. Physical Properties of Heavy and Natural Water 


D,O H,O 
Molecular weight 20-03 18-016 
Specific gravity 1-107 1:0 
Boiling point (°C) 101-4 100-0 
Freezing point (°C) +3-81 0-0 
Maximum density at (°C) 11-2 3-98 
Critical temperature (°C) 371°5 374:2 
Heat of evaporation (kcal/kg) 495-2 533-3 
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Principle of heavy water production 


The production of heavy water is complicated by a 
number of difficulties. Heavy water is present in ordinary 
water, but only in a proportion of | part in 7,000. This low 
concentration means that very efficient methods of separ- 
ation must be used to obtain heavy water of 99-8 per cent 
purity, which is the figure required for nuclear reactors. 

All the production methods in practical use are based on 
an exchange process. It is immaterial whether the enrich- 
ment proceeds directly, i.e. in heavy water itself, or in some 
other deuterium compound. 


Figure 1. Separation factor of deuterium exchange 


If two phases, for instance a liquid and a gas, are 
brought into intimate contact in a container (Figure 1), as 
they may be, for example, by vigorous shaking, and if each 
of these phases consists of a mixture of a hydrogen and the 
same deuterium compound, an exchange or cross-over of 
some of the deuterium atoms will take place. When equi- 
librium is restored, there will be more deuterium in the 
liquid than in the gas. The reason for this—in the case, for 
example, of liquid hydrogen and deuterium—is that, 
owing to its higher vapour pressure, the hydrogen vapo- 
rizes rather more easily than the deuterium. 

It is not essential, however, that the two phases should 
be the vapour and liquid of the same mixture. There are 
other partners which are quite suitable for deuterium 
exchange, as for instance—to name only the most impor- 
tant for technical purposes—water and hydrogen, water 
and hydrogen sulphide, and ammonia and hydrogen. In 
most of these cases, however, exchange takes place only 
in the presence of a catalyst. 

If the proportion of deuterium to hydrogen atoms in the 
liquid is x, and that in the gas y, the exchange process is 
characterized by the ratio « = x/y. Therefore, « is the 
simple separation factor; it indicates the deuterium enrich- 
ment that can be achieved by shaking the mixtures once in 
a container. 
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(1-5 atm) 


Figure 2. Separation factors of various methods of deuterium 
production 


In the diagram shown in Figure 2 the separation factor 
is plotted as a function of temperature for a number of 
exchange processes. The considerable differences in « for 
the various pairs of mixtures is apparent and the marked 
influence of temperature will also be noted. While separ- 
ation factors of 10 are attainable with electrolysis, the 
factor for simple water rectification is only 1-05. 

At first glance, it would naturally seem preferable to 
choose the method offering the highest separation factor. 
The economy of the process, however, is affected by a 
number of other considerations. Thus the power require- 
ments for dissociation of water by exclusively electrolytic 
methods are relatively high. Clusius indicated, as early 
as 1949, a specific power input of 120-150 kW.hr/g of 
heavy water, whereas the plant at Ems which will be 
described below requires only about 2-3 kW.hr/g. 

Other methods have to be carried out either at very high 
pressures or at very low temperatures, or else extremely 
toxic, corrosive, or explosive gases have to be handled— 
all factors which affect the technical practicability of the 
methods. Optimum conditions may be attained by the 
concerted use of several methods: but their attainment 
also depends on whether a heavy water plant can be com- 
bined with other manufacturing processes making use, for 
instance, of hydrogen. At Ems three methods are employed 
in conjunction: electrolysis, rectification of hydrogen and 
deuterium, and rectification of natural and heavy water. 

The extraction of 99:8 per cent heavy water from 
natural water with an initial concentration of only 1/7,000 
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requires an overall separation factor of 3,500,000. The 
figures of 1-05 or 10 attainable in a single step are puny by 
comparison, and must be multiplied many times to enable 
the desired goal to be reached. This is done in rectification 
columns. 

The liquid phase is fed to the top of a column filled with 
packings (Figure 3). An equal number of hydrogen atoms 
in the gaseous phase must now rise in counterflow from 
the bottom of the column. Owing to the large contact 
surfaces created by the packings the equilibrium condition 
defined by the separation factor « is established at every 
point of the column, i.e. at any cross-section the deuterium 
concentration is y in the gas and x in the liquid, x exceeding 
y by the margin of the separation factor. The liquid stream 
thus contains more deuterium atoms than the gas stream, 
and deuterium is therefore carried downwards. Whether 
enrichment takes place or not depends only on what hap- 
pens to the two streams at the bottom of the column. In 
the rectification process the liquid that trickles down is 
quite simply re-evaporated in a vessel known as a still-pot 
or reboiler. For a time, therefore, no product is drawn off. 
Also, since the deuterium is carried downwards without 
anything being removed, it must obviously accumulate at 
the bottom of the column. 

On the right-hand side of Figure 3 the concentration is 
plotted against the height of the column. The diagram 
indicates how the concentration gradually begins to rise 
from the bottom end upwards until the equilibrium state 
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Figure 3. Method of operation of an exchange column 
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represented by the chain-dotted line is reached. The con- 
centrated product can now be withdrawn at the bottom. 
It is most important, however, that no more product is 
withdrawn than is actually being carried down the column. 
If this should happen, the concentration would fall; in 
fact, the concentration at the bottom of the column can be 
controlled by the amount withdrawn. 

These are the essential facts about heavy water produc- 
tion. The rest depends only on accessories, that is to say 
on the equipment required for carrying out the exchange 
process. Surprisingly’ enough, the rectifying columns 
proper form only a very small part of the whole installation. 

Figure 4 shows the three production phases of the plant 
at Ems. The electrolysers of the ammonia synthesis plant 
are connected in cascade for heavy water production. 
Somewhat more water is fed to each stage than is dissoci- 
ated into the components oxygen and hydrogen by the 
electric current. The excess water is condensed out of the 
gases produced and is supplied to the next stage as feed 
water. In this way approximately seven-fold preliminary 
enrichment is attained. The preconcentrate, containing 
about 0-1 per cent D,O, is split up in its entirety into 
oxygen and hydrogen in the last electrolyser. The hydro- 
gen, of which about 400 standard cubic metres is produced 


H,O 
+ 
0-015 % D,O 


per hour, is stored in an intermediate gas-holder and is 
later cooled to about —250°C in the low temperature 
plant. After liquefaction at this temperature, it is con- 
ducted to the low temperature rectifying column. At these 
low concentrations the deuterium atoms occur only in 
combination with hydrogen atoms in the form of hydrogen 
deuteride (HD). Since the equilibrium between H,, HD, 
and D, is ‘frozen’ at this low temperature, the best 
possible product of rectification is 100 per cent hydrogen 
deuteride, which would correspond to 50 per cent Dy. 
In the plant described here the process is carried only as 
far as 60 per cent HD, after which the product taken from 
the column is burnt in a simple oxyhydrogen burner. The 
resulting 30 per cent heavy water preconcentrate is then 
led to a final rectifying column which rectifies water and 
heavy water under vacuum at about 60°C and yields the 
desired product of 99-8 per cent purity. 

This combined process offers the advantage that the two 
rectifying columns can be operated independently of each 
other. The final column can also be used for regenerating 
any highly concentrated heavy water waste without undue 
mixing losses. Another reason for the choice of the com- 
bined process was that the simple and reliable method of 
water rectification was already known and proved, so 
that the risk involved in the first execution of a plant of this 
type was greatly reduced. 


Method of operation of the low temperature plant 
Figure 5 shows the low temperature installation proper 


> H, for the rectification of the liquid hydrogen at —251°C 
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Figure 4. Diagram illustrating heavy water production 99:8 % 
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Figure 5. Diagram of the H,-HD rectifying plant 


(22°K). On the extreme right of the illustration is the 
rectifying column of the Kuhn type. It consists of a series 
of tubes connected in parallel in which the rectification 
process takes place. Refrigerating machinery is needed for 
operating the column. It supplies the refrigeration required 
for cooling the column and for liquefying the hydrogen 
and consists of a compressor, heat exchangers, and 
expansion turbines. 

A three-stage oil-free piston-type compressor brings 
the hydrogen gas up to a pressure of about 14 atm (abs). 
It is then cooled down to liquefaction point by the return- 
ing gas in a series of heat exchangers, each followed by an 
expansion turbine in which it gives up energy to the 
turbine shaft. The last heat exchanger takes the form of an 
ordinary Joule-Thomson exchanger, as used in plants for 
the liquefaction of air. The hydrogen, which is already 
somewhat moist, is now introduced into the bottom part 
of the column to heat the still-pots of the individual 
rectifying tubes and is liquefied while giving up heat to 
these tubes. The liquid is collected and sub-cooled in a 
further heat exchanger before being expanded and intrc- 
duced into a container at the top of the column. From here 
the liquid hydrogen trickles down, evenly distributed by 
capillaries among the various rectifying tubes, and absorbs 
deuterium in contact with the rising vapour. Reaching the 
bottom of the tubes, it is again evaporated by the con- 
densing hydrogen on the outside. It then flows back 
through the tubes in countercurrent to the liquid, is with- 
drawn at the head of the column, and returned through all 
the heat exchangers previously mentioned to the suction 
side of the compressor. Consequently, it first sub-cools 
the liquefied hydrogen and later gives up its refrigeration to 
the fresh incoming hydrogen, absorbing heat in the pro- 
cess, until it is once more at room temperature. This com- 
pletes the refrigeration cycle proper. 
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It will be noted that hardly any refrigeration is given 
up to the outside, with the exception of plant losses. The 
same amount of hydrogen as is liquefied in the heating 
jacket at the bottom of the column evaporates at the same 
time in the inside of the tubes. Although over 1,000 1. of 
liquid hydrogen are produced in the plant every hour, a 
comparatively small refrigeration capacity is sufficient. 

The product, now rich in deuterium, is withdrawn 
through capillaries at the bottom of the tubes. A corre- 
sponding amount of deuterium must be constantly sup- 
plied to the various tubes halfway up the column. This is 
the pre-enriched feed hydrogen coming from the electro- 
lysis plant. As yet, however, it contains many impurities 
and therefore has not only to be cooled but also to undergo 
a purifying process before it enters the column. It is first 
compressed to about 3-7 atm (abs) in a feed compressor, 
after which oxygen is extracted from it in a catalytic con- 
tact furnace. This oxygen then combines with hydrogen to 
form water. The feed hydrogen afterwards goes to two 
interchangeable heat exchangers in which it is cooled by 
returning gas to about — 100°C, so that practically all the 
water is frozen out of it. Only one of these exchangers is in 
operation at any time, while the second is being thawed 
by the fresh incoming gas. After having been cooled to 
about — 200°C in a further exchanger, the feed flows to 
another interchangeable pair of exchangers, as the satu- 
ration temperature of nitrogen has now been reached. 
Electrolysed hydrogen may contain up to two parts of 
nitrogen per thousand, and this is now frozen out in one 
of the last two exchangers. When a sufficient amount of 
nitrogen has been collected, the second exchanger comes 
into operation, while the first is thawed out and the nitro- 
gen led away in the liquid state. In the last stage of purifi- 
cation the whole feed flow is liquefied and re-evaporated. 
The feed liquefier consists, essentially, of a simple tube on 
the outside of which the hydrogen condenses. The liquid 
is then expanded and introduced into the inside of the 
tube, where it again evaporates, leaving the last impurities 
behind. The feed liquefier prevents any obstruction of the 
feed capillaries and the column. 

The overall balance of the process shows that about 
400 m* of hydrogen are introduced through the feed com- 
pressor per hour, of which 1 m%/hr is withdrawn at the 
bottom of the column, while 399 m° flow back through the 
feed exchanger to join the hydrogen used for refrigeration 
and to take part subsequently in the ammonia synthesis 
process. 

The process extracts a small amount of deuterium froma 
flow of hydrogen serving another manufacturing purpose 
and therefore takes ‘parasitic’ advantage of the electro- 
lytic enrichment, the costs of which are debited to the 
ammonia synthesis process. 

At the bottom of the diagram, purifying equipment is 
shown; this is required for the first commissioning of the 
installation. It is of the greatest importance that the gas 
flowing in the refrigerating circuit should be absolutely 
pure, as otherwise the expansion turbines might become 
choked with ice. When the installation receives its first 
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Figure 6. Temperature-entropy diagram of the refrigerating process 


charge, the hydrogen is therefore cooled with liquid air 
and cleaned in an activated carbon filter. A dry gas-holder 
stores up this pure charge when the plant is out of opera- 
tion. 

The refrigeration cycle illustrated in Figure 5 is plotted 
in the temperature—entropy diagram in Figure 6. The first 
expansion of the gas after three-stage compression takes 
place in two turbines connected in series, so that in all 
there are four turbines, two of which operate below the 
temperature of liquid air and one only about 30°C above 
absolute zero. The installation differs from known gas 
liquefying or separating systems in that all the turbines are 
connected in series, so that the total amount of gas handled 
by the circulating compressor flows through each of them. 

We will next discuss the component parts of the heavy 
water separation plant. 


Compressors 


Both compressors operate without any oil contamina- 
tion.! The piston rod is guided at two points outside the 
cylinder, in which the piston moves freely, sealing being 
obtained by labyrinth grooves. The stuffing-box is also 
provided with labyrinth rings, so that the gas handled 
cannot come into contact with oil at any point. As no 
stuffing-box losses are admissible, the whole crankcase of 
the compressors at Ems is made gas-tight and is connected 
to the gas-holder by way of an oil trap. In practice this 
trap has hardly any work to do, since the pressure in the 
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gas-holder is constant and the gas in the crankcase there- 
fore remains unchanged. An intermediate tapping at the 
stuffing-box ensures that leakage losses are returned direct 
to the suction side of the cylinders and do not enter the 
crankcase at all. 

Figure 7 (left) shows the feed compressor. There is a 
comparatively large valve station (not shown in the Figure) 
on the suction side of the machine, its size being explained 
by the fact that it is necessary to aspirate either nitrogen 
or normal or enriched hydrogen, according to whether 
the plant is in normal operation or is being flushed with 
nitrogen for overhaul purposes. 

The circulating compressor of the refrigerating plant 
(Figure 7 (right)) is a three-stage machine, as already 
mentioned, and is directly coupled to a 3,000 V driving 
motor. This compressor also has a completely enclosed 
crankcase. 

The use of oil-free compressors has considerably simpli- 
fied the design of the low temperature plant, since no 
equipment is now required for removing oil mist down- 
stream of the compressors and there is no danger of the 
heat exchangers and the column being fouled. 


Vacuum insulation of the heat exchangers 


One of the main problems raised by the design of the 
plant was the layout and insulation of the heat exchangers 
and control valves. It was not possible to use mineral wool 
as in ordinary air liquefaction plants, since in two-thirds 
of the installation the prevailing temperatures are below 
that of liquid air. However thick the layers of insulation 
were made, the air contained in them would condense on 
the inner walls, and there might even be a dangerous 
build-up of oxygen, involving the risk of explosions if the 
slightest hydrogen losses should occur. On the other hand 
it is uneconomical to fill the insulation space with hydro- 
gen, as the thermal conductivity of this gas is about seven 
times that of air. 

To circumvent all these difficulties, vacuum insulation 
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was chosen. The designers were aware that it would not be 
a simple matter to make such large pieces of equipment 
tight against high vacuum, especially if they contained 
hydrogen under pressure. A vacuum of at least 10~* torr, 
or a pressure of about 1/10* atm, is needed for good insu- 
lation. On the other hand, the advantages were so attrac- 
tive that the attempt seemed justified. With high vacuum 
insulation, when the mean free path of the molecules is 
longer than the distances between individual units of 
equipment these distances cease to be important, and the 
whole plant can be made very compact. In addition, it is 
not necessary to cool tons of insulating material when the 
temperatures are being lowered prior to operation, and 
starting times are therefore much reduced. 

A second basic decision affecting the design of the plant 
was that all regulating valves should be remote controlled 
with gas under pressure. It has hitherto been the custom in 
low temperature installations to connect the cold valves 
to hand-wheels outside the outer jacket by means of 
spindles of suitable length. This severely restricts the free 
placing of the control valves, while the numerous spindles 
cause additional heat and leakage losses. With servo- 
controlled valves, whose servomotors are also at low 
temperature, only two small calibre pipes have to be 
carried to the outside, while the valve can be placed exactly 
as desired. 

Figure 8 shows a cross-section of the vacuum jacket 
containing the column. There are four such units in the 
whole plant. The first contains the heat exchangers down 
to — 100°C, the second the exchangers down to — 200°C, 
the third the feed liquefier and final feed heat exchangers, 
and the fourth the column. 

All the vacuum jackets are of the same fundamental 
design. The heat exchangers and the column are all sus- 
pended from long thin stainless-steel rods in a lower space. 
As stainless steel has a very poor thermal conductivity, 
undesirable heat gains can be kept to a minimum. The 
jacket is connected to an oil diffusion pump through a 
lateral branch. There are no flange connections or other 
mechanically detachable joints in the jacket itself. All 
pipes are welded or soldered. All components incorporat- 
ing flange connections, and particularly the control valves, 
are mounted under a separate hood in the upper part of 
the unit. In this upper space a rough vacuum only is main- 
tained, since there is obviously a possibility of hydrogen 
leaks developing in some of the flanges in the course of 
time. To ensure good insulation, however, this upper hood 
is of double-walled design, and the space between the 
walls is kept under high vacuum as in a Dewar vessel. 

The turbines are located laterally on the suction pipe 
of the diffusion pumps and are also covered by a vacuum 
hood. Radiation shields of aluminium are fitted inside all 
the vacuum jackets in order to reduce radiation losses 
to a minimum. For the purpose of erection and dismant- 
ling, the large vacuum jackets can be simply lowered and 
the hoods over the fittings removed from above. The 
various units of equipment are then readily accessible 
from the platforms. The lowering of the vacuum jackets 
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takes much less time than would the removal, for instance, 
of slag-wool insulation. 


Servo-controlled valves for ultra-low temperatures 


Extensive development work made it possible to design 
servocontrolled changeover and regulating valves which 
operate reliably even at extremely low temperatures 
(Figures 9-11). The sectional drawing in Figure 9 shows a 
three-way valve with low temperature servomotor. The 


Figure 8. Section through a vacuum jacket with column, heat ex- 
changers, diffusion pump, regulating valves, and expansion turbine 
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Figure 9. Three-way valve with low temperature servomotor 


operative component of the servomotor is a bellows unit 
that transfers its force through a bell to the central spindle. 
This in its turn actuates the valve disk. All seals are 
metallic and are tightened with extensible shank bolts 
which ensure adequate sealing pressure even when big 
temperature changes occur. The servomotor spring can be 
fitted to operate in either direction, so that the valve always 
takes up the desired position when free of pressure. Hydro- 
gen is of course the only choice for the actuation of the 


Figure 10. Regulating valve 
for liquid hydrogen 
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servomotor, as air would solidify. Even with hydrogen 
there is some risk of its liquefying when the servomotor 
operates at a very low temperature level. For this reason 
the servomotor is separated from the valve by a thin walled 
stainless-steel cylinder and is warmed by a gas heating 
system enough to ensure that the control gas cannot 
condense, although no great losses can occur. 

The regulating valves (Figure 10) are similarly designed. 
The one shown in the Figure is for liquid hydrogen. 
Whereas in the three-way valve two end positions only are 
needed, a regulating valve must take up any desired 
position, and this must be exactly adjustable from outside. 
To this end the whole valve is made absolutely frictionless. 
The spindles are held by elastic springs, and the infinitely 
adjustable control pressure acting on the bellows operates 
against a spring, so that the pressure directly determines the 
spring loading and thus the valve position. Here again the 
servomotor is divided from the valve by a distance piece 
and separately heated. 


Expansion turbines for hydrogen 


The most important units in the refrigerating plant are 
the expansion turbines. Without these no refrigeration 


Figure 11. Valve assembly from the side. The cold end of the vacuum 

jacket carries the valves and is itself supported on the warm tank 

carrier by means of a thin walled stainless-steel cylinder. The heating 

coil visible in the Figure is used for heating the whole valve space when 

overhauls are carried out. The control and heating pipes for all servo- 
motors are also shown 
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would be produced, and all the other parts of the plant 
thus depend on their efficient operation. 

The installation at Ems is the first hydrogen liquefying 
plant in which refrigeration is produced exclusively by 
expansion turbines. Most of the refrigeration required for 
hydrogen liquefaction is normally obtained from the 
evaporation of liquid nitrogen, which means that the 
hydrogen liquefier has to incorporate a complete liquid 
nitrogen plant. 

The turbines (see Figure 12) are, like most expansion 
turbines for low temperature applications, of the centri- 
petal type, i.e. the gas flows radially inwards through the 
turbine wheel and leaves in the axial direction. The turbine 
casing is supported on the warm fixing flange by means of 
a thin stainless-steel cone, and an internal bearing housing 
encloses the whole oil-lubricated part of the turbine rotor. 
The only insulation between the cold turbine casing and 
the warm bearing housing is provided by high vacuum. 
The turbine rotor, with a vertical shaft, runs in two oil- 
lubricated bearings and has an oil brake in the centre to 
dissipate the excess energy. The turbine wheel is overhung 
and projects into the cold casing. A magnetic pick-up 
excited by the bottom end of the shaft measures the turbine 
speed. 

When the turbines are fitted, the gas entering them is 
first passed through a filter, in order to prevent damage to 
the inlet and turbine wheels. A hand valve is provided on 
the inlet and outlet sides to permit the turbines to be 
removed singly even when the rest of the plant is cold and 
contains hydrogen. 

The design of expansion turbines for hydrogen presents 
several difficulties. First of all a seal must be provided for 
a shaft turning at about 100,000 rev/min. The distance 
between the cold turbine wheel and the warm lubricated 
bearing is only about 30 mm, yet the temperature differ- 
ence is over 250°C. It is imperative that no trace of oil 
whatsoever should be able to find its way into the turbine, 
and refrigeration losses through the casing and shaft 
must also be kept low. 


Figure 12. Cross-section through the expansion turbine for hydrogen 
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Figure 13. Characteristic of the expansion turbines for hydrogen 


A further difficulty arises from the general turbine 
characteristic (Figure 13). The efficiency reaches its 
optimum at a figure of about 065 for the ratio of the 
peripheral speed of the turbine wheel U to the head 
velocity Cy. This is the velocity attained when the gas is 
isentropically expanded with the use of the full turbine 
head. Even at a pressure ratio of about 2, C, reaches the 
speed of sound, which for hydrogen at room temperature 
is over 1,200 m/sec. The corresponding optimum periph- 
eral speed would thus be 800 m/sec, a figure which cannot 
be reached with the materials so far available. This means 
that only a comparatively low head can be attained in one 
turbine stage, even if peripheral speeds up to 400 m/sec are 
adopted, and this is one of the reasons why the plant is 
equipped with four turbines, whereas normal air lique- 
faction plants have one turbine only. 


Hydrogen-hydrogen deuteride rectifying column 

Figure 14 shows the low temperature installation with 
the vacuum jackets lowered. The rectifying column for 
liquid hydrogen, consisting of 90 parallel tubes, hangs in 
the foreground. About two-thirds of the length of the 
column is used for the rectification process proper, while 
in the bottom third the hydrogen evaporates as it trickles 
down through the tubes. The hydrogen condenser sur- 
rounding the tubes also serves as a liquid container. The 
feed hydrogen is supplied to the individual tubes through 
capillaries about half-way up. At the head of the column, 
thin feed tubes can be seen projecting into each tube of the 
column. These contain the capillaries for the reflux liquid. 
The upper liquid tank is hidden under the hood. It has an 
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extension in the centre to take the level measuring equip- 
ment. Reflux in the column is kept constant by exact 
control of the liquid level in this tank. The last heat 
exchangers for partial condensation of the circulating gas 
are also visible above the column. 


Vacuum plant 


The vacuum plant (Figure 15) embraces three different 
systems. The system generating the high vacuum in the 
main vacuum jackets comprises a mechanical rough 
vacuum pump and aRoots blower, both connected to a 
header to which the four oil diffusion pumps are joined up 
in parallel. The high vacuum is therefore generated in 
three stages. A second system again includes a mechanical 
rough vacuum pump with a single diffusion pump con- 
nected to it and operating on a high vacuum header. The 
spaces between the double walls of the valve hoods and the 
hoods over the turbines are connected to this header. A 
third system with a single rough vacuum pump evacuates 
the inside of the valve hoods, where a less stringent 
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Figure 14. Rectifying column for liquid hydrogen 
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Figure 15. Diagram of the vacuum plant 


vacuum is required. The three normal pumps are supple- 
mented by a stand-by which can take over for any of the 
other units through a header, and can also pump out any 
vacuum jacket that has been shut down by way of a 
separate pipe before it is connected up to the main header 
again. Each of the systems has its own vacuum monitor 
which isolates all the jackets by pneumatic means as soon 
as the pressure rises above a certain value. At the same 
time all the electrical equipment, such as the heating system 
of the diffusion pumps and the electric vacuum measuring 
units, is switched off to prevent any oxyhydrogen gas 
present in the vacuum space from being ignited by the 
electric current. 


Control room 

The whole heavy water plant is supervised from a central 
control room. A diagram of the plant, as shown in 
Figure 5, appears on one of the walls. Hand-wheels are 
fitted to correspond with the positions of the valves in the 
diagram, which greatly simplifies control of the plant. 
Automatic recording apparatus for temperatures, pres- 
sures, liquid levels, and vacua can be seen in the back- 
ground, and on the left the revolution indicators for the 
turbines. 


Hydrogen deuteride combustion 

The combustion equipment takes up less than | m? of 
floor area. It comprises a burner followed by a flame tube 
of high temperature steel, a cooler, several changeover 
valves, and a cold trap for freezing the heavy water out of 
the outgoing excess oxygen. Here again a safety system 
ensures that no oxyhydrogen explosions can occur. 


Heavy water rectification 


Figure 16 shows the heavy water rectification column. 
The liquid flowing down from the burner is collected in a 
container and fed by way of a float regulator to the feed 
evaporator. The feed is metered by setting the desired 
heating performance, and the feed vapour is then distri- 
buted through capillaries to the individual rectifying 
tubes. Here, heating of the still can be done with ordinary 
steam, and the reflux is produced at the head of the column 
in a simple water-cooled condenser. From here it goes to a 
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vessel whose overflow level determines the reflux quantity. 
The overflow of this vessel is the head product. It is col- 
lected in a separate receiver and returned to the electro- 
lysers. The heavy water is withdrawn in the form of steam 
through capillaries at the bottom of the tubes. Removal is 
automatically controlled by a concentration regulator. 
This employs a simple glass float which is calibrated to the 
specific gravity of heavy water and is scanned by photo- 
cells. When the concentration begins to fall, the float moves 
downwards, causing a reduction in the amount of product 
withdrawn. About | 1./hr of 30 per cent heavy water in all 
is introduced into the final column per hour, and is there 
separated into the ‘bottom’ product of 99-8 per cent 


i 03U./hr 


1. Reflux condenser 2. Kuhn column 

3. Feed distributor 4. 30% D.O feed 

5. Feed evaporator 6. Vapour 

7. Insulation 8. Heating jacket 

9. Product condenser 10. Concentration regulator 
11. Product (99-8% D,O) 12. Manostats 
13. Vacuum pump 14. Head product (3% D,O) 


Figure 16. Diagram of the heavy water rectification system 
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purity and a ‘head’ product of 3 per cent concentration. 
The vacuum in the column and in the product condenser 
is automatically controlled by manostats. 


Conclusion 


A few comments on the time required for the develop- 
ment of this plant may be of interest. It was necessary first 
of all to test all the component parts singly, a procedure 
which was in itself fairly time-consuming. 

Wide use was made at this time of the hydrogen liquefier 
in the Physico-Chemical Institute of Ziirich University, 
which was needed for cooling down the low temperature 
fittings, testing soldered joints, and calibrating thermo- 
couples. In the same Institute the efficacy of the packings 
used for separating hydrogen and deuterium was first 
tested and the feasibility of rectifying hydrogen with Kuhn 
columns was proved. 

Before the final decision was taken to insulate the whole 
plant with high vacuum, a single vacuum-jacket unit was 
erected in the open at Oberwinterthur. The same unit, 
with the same heat exchangers, was later erected practi- 
cally unchanged in the installation at Ems. Immediately 
after assembly a vacuum of 10~‘ torr was attained, so that 
the decision to carry out the plant with high vacuum 
insulation at once appeared justified. At Ems it was later 
possible to generate vacua of up to 10~® torr in the big 
jackets. 

The experimental plant was cooled down to the tem- 
perature of liquid nitrogen with a first prototype turbine, 
so that it was now clear that all the components could be 
relied upon even at very low temperatures. 

It may be mentioned here that the technical realization 
of the whole plant demanded great persoral efforts from 
all concerned, since hardly a dozen persons were available 
for the calculation and designing work throughout the 
development period of just under four years. 


We should like in conclusion to express our special 
thanks to Professor Clusius of Ziirich University, to Pro- 
fessor Grassmann of the Swiss Federal Institute of Tech- 
nology in Ziirich, and to Professor Kuhn of Basle Uni- 
versity for their valuable assistance, as well as to the 
Management of the Emser Werke for their.close and 
helpful co-operation. 
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THE electronic specific heat of a superconductor in the 
normal state may be calculated from the observed critical 
field curve. Similarly, the change in electronic specific heat 
with pressure may be calculated from the change in critical 
field under pressure. Such work on a wide range of super- 
conductors! has shown that the electronic specific heat is 
much more sensitive to volume changes than the Sommer- 
feld free electron theory predicts. 

It would obviously be interesting to discover if this 
volume sensitivity can be understood from a study of the 
actual electronic structure of a metal. To provide.a basis 
for such an investigation we have thought it useful to col- 
lect experimental information on the pressure dependence 
of the critical field in aluminium. This is the supercon- 
ductor about whose Fermi surface most is known. The 
de Haas—van Alphen effect in this metal has been studied 
extensively by Gunnerson,? and work on cyclotron reso- 
nance and on the anomalous skin effect is also available. 
Detailed theoretical studies have been made by Heine* 
and by Harrison.*® 

Grenier® and Muench’ have measured the pressure 
dependence of the transition temperature 7, in aluminium. 
A calculation of the pressure dependence of the electronic 
specific heat requires a knowledge of the effect of pressure 
on the critical field at lower temperatures as well, and we 
have now made such measurements at temperatures down 
to around 0-3°K. It is found that the electronic specific 
heat is highly volume sensitive in this metal too. 

In addition to this work on aluminium we have also 
examined the effect of pressure on the critical field of zinc. 
Here the technique used for creating the pressure does not 
yield a pressure which is sufficiently isotropic to allow 
quantitative conclusions. We include the results here for 
completeness. 


Experimental method 


Adiabatic demagnetization of ferric ammonium alum 
was used to cool the specimens. The salt was mixed witha 
cold-setting resin (‘Prestolith’, obtainable from Althaus 
and Vogt GmbH, Bochum-Harpen) and compressed 
under a pressure of 500 kgcm~? around copper fins. This 
produced a remarkably stable pill which has now been in 
use for about 18 months without apparent deterioration. 
According to Andres® the thermal resistance between 
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copper and salt is of the order of 2 x 10* deg.cm? W~! at 
0-4°K. 

The pressure was made by the ice-bomb technique first 
used by Lazarev and Kan? in similar experiments on tin 
and subsequently by several other investigators. The 
specimens to be put under pressure are enclosed in a thick 
walled container which is filled with water and which can 
be firmly closed. On cooling to liquid helium temperature, 
pressures of about 1,700 atm are developed. 


Specimens 


X-Y recorder 


Amplifier 


Figure 1. The apparatus 


The ice-bomb contained three specimens: one each of 
aluminium and zinc, and also one of another metal used 
for pressure calibration. Three other similar specimens 
were mounted in a similar container without ice, and both 
containers were connected thermally to the copper vanes 
in the salt pill. Both containers were therefore at the tem- 
perature of the salt. The temperature adjustment in the 
pressure container was excellent; in the container without 
ice it could only be maintained by soldering the specimens 
to the container. 

The critical fields were measured in the following man- 
ner. A magnetic field which increased linearly with time 
was applied, and the difference in the rate of change of 
flux in two coils wound round the containers was observed. 
This was done by amplifying and recording the output 
from the coils which were wound in opposition. Figure | 
shows a schematic diagram of the apparatus, and Figure 2 
shows typical recordings. When nothing happens, the 
output voltage is zero, but deflections occur when one or 
the other of the specimens becomes normal. The rates of 
increase of field could be varied between 0-1 and 0-5 gauss 
sec~!. Too rapid an increase led to smeared-out transitions, 
and usually a rate of change of 0-15 gauss sec~! was used. 
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Figure 2. Typical recorder traces; the downward deflections come 
from the specimen at zero pressure 


The first measurements were made using elliptic speci- 
mens 3 mm in diameter and 30 mm long, but it was found 
that the transitions were too wide to be properly inter- 
preted. The measurements shown here were therefore 
made with specimens in the form of wires | mm in diameter 
and 30 mm long. The very small indium joints used to 
make thermal contact between specimens and container 
caused a local distortion of the magnetic field. This was 
too limited in extent to affect the results. 

The pressures were measured in each run by means of a 
specimen with known 0H,/dp (H, is the critical magnetic 
field, and p is the pressure). Alekseevskii and Gaidukov’® 
have made measurements on cadmium, and we used this 
substance as a pressure calibrator in our first experiments. 
Unfortunately, it now appears that the pressure effect in 
cadmium is sometimes sufficiently anisotropic to give 
highly misleading results for the pressure in the ice-bomb. 
This led to an overestimate of 0H,/dp at low temperatures 
in preliminary measurements on aluminium reported 
briefly elsewhere." In the work reported here we used 
indium to measure the pressure. Indium is also anisotropic, 
but this metal appeared to give consistent values of the 
ice-bomb pressure. This is presumably because indium 
flows so easily under pressure that anisotropic stresses 
cannot build up. Muench’ also confirms the reliability of 
indium for ice-bomb calibration. 


Results 


The critical field H(p) at pressure p is given approxi- 
mately by 


A Ap) = Hop) t1-(T/TAp)P} ...() 


where H,(p) is the critical field at the absolute zero, and 
T(p) the critical temperature at pressure p. The plot of 
H(p) against T? should therefore be approximately a 
straight line. Similarly a plot of AH, =H (0)—H Ap) 
against 7? should yield a straight line. 

Figure 3 shows a graph of H.(p) against T? for alu- 
minium with p = Oand p = 1,700 + 100 atm. There appears 
to be some deviation from the linear behaviour predicted 
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by equation (1). This happens in many superconductors, 
but is surprising in this case since Goodman and Men- 
doza!2 found that equation (1) is obeyed to better than 
1 per cent for p = 0. Our temperature determination may 
be less accurate than theirs, and this point requires further 
investigation. 


(gauss) 


Critical field -——»> 


0 05 0 rade 
Tobe (°K)? 


Op=0 @ p=1,700+ 100 atm 
Figure 3. The critical field of aluminium 


Figure 4 shows AH, for two pairs of aluminium 
specimens and for two pairs of zinc specimens. Three to 
four demagnetizations were made with each pair of 
specimens. The relative accuracy in AH, is lower than in 
H., and possible deviations from a linear dependence of 
AH, on T? cannot be seen. The pressures in the ice-bomb 
were approximately equal in all runs, i.e. p = 1,700 + 100 
in all cases. 


(gauss) 


4h.=—> 


0 0:5 1-0 S 
¢— (°K)? 


AVY Zinc 
Figure 4. AH. in aluminium and zinc; p =1,700+100 atm in all cases 


© Aluminium 


One pair of aluminium specimens was made from 
spectroscopically pure material obtained from Johnson 
Matthey and Co. Ltd. The other pair was made from zone- 
refined material. (We are grateful to Dr. W. Rohner of the 
Research Laboratory of the Aluminium Industrie A.G., 
Neuhausen, for making this available to us.) In both cases 
the specimens were drawn to size and left unannealed. In 
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aluminium the same values of AH, were obtained in all 
runs, and a definite value for 0H,/dp can be given. We find 


OH,/dp = —[3-1+0-2+ 
, + (1:7 + 0:2) 27] x 10-® gauss dyn-! cm? 


where t = 7/T,. 

All the zinc specimens were cast from spectroscopically 
pure material from Johnson Matthey and Co. Although 
the results on AH, obtained for each pair of specimens are 
consistent, and can be plotted as a linear function of T?, 

different specimen pairs give different values of AH,. We 
believe that this must be caused by the presence of aniso- 
tropic stress systems within the ice-bomb, and that each 
set of results corresponds to a different stress system. If the 
width of the transition can be used as a guide then the 

results on curve A may perhaps be thought to correspond 
to the most nearly hydrostatic pressures. This curve yields 
for zinc 


0H,/op = —[1-4+0-34+ 
+ (0-4 + 0-2) ¢?] x 10-® gauss dyn-! cm? 
Curve B gives 


0H,/Op = —[0-6+0-24 
+(0:5 + 0-2) 77] x 10-4 gauss dyn-! cm? 


On one occasion, badly smeared-out positive values for 

_0H,/Op in zinc were obtained. 

_ In cadmium, our results agreed with those of Alek- 
seevskii and Gaidukov”® in some runs. In others we found 
a change in critical field three times smaller than theirs. 

Presumably, then, 0H./dp depends on the orientation of 
the stress in this metal too. 

It is obviously possible to calculate the pressure de- 
pendence of 7, from these measurements, and we have 
pointed out elsewhere? that dy/dp may be calculated, too. 
(yT is the electronic specific heat per mole.) If H, = H,.f(t), 
where ¢ = (7/T,), and if the form of f(t) is independent of 
pressure, then 

Wee = af) — yf (lt 
Pp op 


1 ay a “n 


+4H)tf'(t) |: ee a 
where f(t) is the derivative of f(t) with respect to t, and v 
is the molar volume. This may be expressed somewhat 
more usefully in the form 


diny _ _,[(4./@P) m0 ¢, | Gut 
ca [a ees ding 


We denote the quantity dIny/dInv by g in the following. 
Clearly also 


dinT, K (0H, OH \= 
pay Hea HOR 
dinw t Op } r= Te OT | Le 


where K is the bulk modulus. 
In calculating g it is found that the results are sensitive 
to the value of f’(1) used in equation (3). In the case of 
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aluminium, the measurements by Goodman and Mendoza 
indicate f’(1)= —2, while our work suggests that 
f'(1) = — 1-85 may be possible. A final decision must await 
a better determination of the critical field curve of alu- 
minium. The results of such calculations are collected in 
Table 1 together with values of dIn7./dInv. For zinc, 
both curves A and B shown in Figure 4 have been used. 


Table 1 
| 
din 7. diny =g 
dinv : dinv — 
Aluminium: | | 
Muench | 12+1 — 
Grenier | 13+1 — 
This work: f’(1) = —2:0 ent 10+3 
This work: f/(1) = —1-85 a3 
Zinc: | | 
Curve A 6+1 743 
Curve B 4+1 | 1+3 


The probable errors are estimated from the scatter in 
our points. This does not exclude the possibility of 
systematic errors, and it appears! that higher values of g 
are obtained from direct measurements of 0H,/dp than 
from measurements of the volume change on destruction 
of superconductivity. White! has recently used the 
thermal expansion coefficient at low temperatures to de- 
termine g in aluminium. He finds a smaller value for g 
than either of those suggested here. (We are very grateful 
to Dr. G. K. White for informing us of the results of 
preliminary measurements.) 

The measurements are of interest from two points of 
view. The volume dependence of y gives information on 
the behaviour of the density of states N(O) per mole per 
unit range of energy at the Fermi surface in the normal 
metal since 


yT = 3.7n7k? N(O)T 5 tH) 


where k is Boltzmann’s constant. The volume dependence 
of T,, on the other hand, is closely connected with the 
behaviour of the interaction which causes superconduct- 
ivity. We will discuss these points separately. 


The density of states 
For a system of free electrons, N(0) is given by 
2mmv (=) 


NO) = 3 


. (6) 


where m is the electronic mass, n the number of electrons 
per cubic centimetre, and wv the molar volume of the metal. 
This clearly leads to 


g = din M(0)/dInv = 3 263 (7) 
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The values found for g in the present work, just as those 
found for other metals, are much larger than 3. This is 
perhaps not surprising in view of the complicated structure 
of most of the metals investigated. Yet it must be noticed 
that din N(0)/dInw differs much more from the free 
electron value than (0) itself does. 

This problem has recently been discussed by Rohrer.’® 
He showed that small regions of overlap of the Fermi 
surface from one Brillouin zone into another give contri- 
butions to N(O) which are very sensitive to changes in the 
energy gap at the zone boundary. This energy gap itself 
depends upon the interatomic distance. 

In the case of aluminium, Harrison® has pointed out 
that the free electron model leads to a full first zone, a 
second zone almost full of holes, and a third zone with a 
somewhat complicated region of occupied states. In the 
fourth zone, too, there are minute occupied regions. The 
first zone, therefore, contributes nothing to M(0). The 
second zone roughly 80 per cent and the third 20 per cent 
to the total value of N(O), while the pockets in the fourth 
zone might contribute about | per cent. 

More detailed calculations by Harrison,® taking the 
effect of the periodic potential better into account, confirm 
that the first zone is full. (This disagrees with earlier views 
due to Heine.*) It appears, however, that the pockets in the 
fourth zone disappear. The area of the second zone surface 
is almost unaffected, while that in the third zone appears 
to be increased by the improved approximation. The area 
in this third zone is particularly sensitive to changes in the 
periodic potential since much of the surface lies very close 
to zone boundaries. It must be this surface which con- 
tributes to our large value of g. 

In a quadrivalent metal with the same crystal structure, 
the overlap into the third zone is much bigger, and the area 
of Fermi surface in the third zone should be much less 
sensitive to changes in periodic potential. Thus we would 
expect a smaller contribution to g from the third zone ina 
quadrivalent metal than in a trivalent one (for a given row 
in the periodic table). Therefore, unless large contributions 
are made by the fourth zone or unless this argument is 
invalidated by the large difference in atomic number we 
might expect that in a substance such as lead g should be 
smaller than in aluminium. Unfortunately, the available 
experimental data is not yet sufficiently reliable to test 
conclusions of this sort. For aluminium the present work 
gives g x 8. Until recently the corresponding value for lead 
was g = 2, thus supporting the discussion above. Now, 
however, very careful work on lead by Mapother and 
Garfinkel!® has suggested g ~ 6 in lead. 

Our results on zinc show that in this metal, too, the 
density of states is strongly affected by a change in volume. 
It is interesting to note that recent work by Dmitrenko, 
Verkin, and Lazarev!’ on the de Haas—van Alphen effect 
under pressure shows remarkable changes in the period of 
the observed oscillations. This corresponds to a rapid 
change in size of the pockets of charges giving rise to these 
oscillations. These pockets may thus be important in 
contributing to the observed high value of g. 
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The transition temperature 


Bardeen, Cooper, and Schrieffer'® give an expression for 
the transition temperature 7, of a superconductor. This 
may be written in the form 


T. = 0:850yexp[—1/M(0) V] oe 


where Op is the Debye temperature, N(0) is the density of 
states at the Fermi surface at T = 0, and V is an interaction 
constant which determines the strength of the interaction 
causing superconductivity. On examining the result of a 
volume change, we find 


din(T,/9p)/dInv _ din(N(0) V) 


n (0-8. _ ding 


’p)/C (9 
In(0-850p/T.) ed 


It is obviously of interest to investigate the effect of 
volume change on (0) V, and this has been done theoreti- 
cally by Pines!® and by Morel.?° It might be thought that 
this would be closely correlated with the volume depen- 
dence of (0) itself, and Morel’s work tends to support 
this conclusion. Pines, on the other hand, obtains an 
expression for N(O) V which depends only on the valency 
and interatomic distance in the metal. He suggests that 
changes in N(O) should be compensated by changes in V. 
In this case dInN(0)V/dInw should not be correlated 
with g. This is supported by an observation of Rohrer!® 
who finds that 

din N(0)V/dInw ~ 2 


for most soft superconductors. The mutual compensation 
of changes in (0) and in V is especially clearly seen in a 
plot against g of all observations at present available of 
din N(0) V/dInw in soft superconductors. This is shown 
in Figure 5. 
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Figure 5.d\n N(O)V/d\Inv and its dependence on g in the soft 
superconductors ; 
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THE merits of producing deuterium by the distillation 
of hydrogen are quite well known. However, the best 
methods of distillation column design are not so well 
established. 

The usual correlations available for predicting column 
design parameters fail when applied to the hydrogen 
systems, since the bulk of the data for the correlations is 
obtained from higher boiling systems such as aqueous 
solutions and hydrocarbons. For example, the optimum 
downcomer area was found to be much larger than pre- 
dicted, and one well known correlation! gives an overall 
plate efficiency of 130 per cent. 

These and other correlations are generally applicable 
only to systems having physical characteristics similar to 
those for which the data were obtained. The unconven- 
tional properties of liquid hydrogen, which has a density 
one-fourteenth, and a viscosity one-twentieth that of boil- 
ing water, give rise to unusual behaviour in a distillation 
column. Accordingly, research has been conducted at the 
National Bureau of Standards Cryogenic Engineering 
Laboratory in Boulder, Colorado, to determine the design 
parameters for a hydrogen distillation column. 
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Since the behaviour of packed columns?-> and bubble 
cap columns‘ has been previously reported, this research 
was primarily concerned with the characteristics of per- 
forated or sieve plate columns. Such parameters as plate 
spacing, plate geometry, maximum allowable vapour velo- 
city, weir height, and downcomer area were investigated, 
and the results of these investigations were applied to the 
design of a pilot plant column. 

Since deuterium occurs naturally in the form of HD 
rather than D,, the pilot plant was designed to separate the 
system H,-HD. However, all the experimental work lead- 
ing to the design of the pilot plant was carried out with the 
H,-D, system, for D, is readily available, and HD is not. 
The use of the H,-D, system in the experimental work 
rather than the H,-HD system is easily justified. Since the 
two pairs have very similar physical properties, they should 
have similar distillation behaviour. As an index for com- 
parison, the Drickamer and Bradford! correlation predicts 
that the two systems will have the same overall plate 
efficiency. The O’Connell’ correlation predicts a plate 
efficiency for the HD system only slightly higher than that 
for the D, system. 


Figure 1. A schematic flow diagram of pilot plant 
used for low temperature distillation 


pump 


O 
HE. Flowmeters 
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The initial distillation work was done with a series of 
1-Sin. O.D. columns (1-35 in. I.D.), and the results of these 
preliminary experimental studies have already been pub- 
lished.® These studies culminated in the design and opera- 
tion of a hydrogen distillation pilot plant. The design of 
this plant was covered in the same article,s and some 
details of its operation have also been reported.® However, 
the observed plate efficiencies have not been satisfactorily 
explained. This is the subject of the present paper. 


The pilot plant column 

The pilot plant column contains thirty plates 6 in. in 
diameter spaced 4 in. apart, and made of 40 mesh brass 
screen. Maximum feed rate is 14-2 ft3(n.t.p.)/min. The 
column was operated with vapour velocities ranging from 
0-2 to 2-0 in./sec, and reflux ratios from total to an (L/D) 
of 2:05. A schematic diagram of the column is shown in 
Figure 1. A set of data typical of a partial reflux run is 
recorded in Table 1. For details of the design and operation 
consult the previous articles.*® 


Experimental results 

The column has been operated repeatedly with hydrogen 
for periods of 1-5 days of continuous operation. The most 
interesting information obtained from these runs concerns 
the plate efficiency. 


Overall plate efficiency 

The principal object of experimental runs was to 
determine the overall plate efficiency. Sample taps, located 
at seven points along the column (Figure 1), fed gas phase 
samples to a continuously recording thermal conductivity 
analyser.!° These analyses were recorded along with the 
other process variables so that the overall plate efficiency 
could be determined as a function of vapour velocity. This 
dependence is shown in Figure 2. 

These plate efficiencies represent data for runs during 
which the reflux ratio was varied from total to an (L/D) 


Plate efficiency (per cent) 


Superficial vapour velocity (in./sec) 


() 1:35 in. columns, sieve plates A 6:00 in. column, sieve plate 
xX 1:06 in. column, bubble cap (Sellers and Augood) 


Figure 2. Plate efficiencies of three tvpes of column 
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ratio of 2:05, and during which the composition of the 
mixtures on the plates varied from pure hydrogen to 98 
per cent deuterium in hydrogen. No dependence of overall 
efficiency on either reflux ratio or composition could be 
found. 

Three points at low vapour velocity, from 0-1 to 0-5 
in./sec, group themselves around a low value of plate 
efficiency. This behaviour at low vapour velocities is not 
uncommon, and has been reported by several authors!) 
who investigated higher boiling systems. The effect is one 
of hydraulics on the plate, and not one due to peculiarities 
of a given system. It is therefore reasonable to expect to 
find this behaviour even in the hydrogen system. 

Two sets of data, in addition to the pilot plant efficiencies, 
are shown in Figure 2. One set is that of Sellers and 
Augood,® which was obtained with a single bubble cap in 
a 27 mm diameter column. This set correlates very well 
with that of the pilot plant. 

The other data of Figure 2 is that obtained from the 
1:35 in. diameter columns. These plate efficiencies are 
significantly higher than those of the pilot plant and of 
Sellers and Augood. This curious discrepancy will now be 
examined. 


Plate efficiency relations 


On a theoretical plate, by definition,!* 14 the average 
composition of the vapour leaving is the equilibrium value 
for the liquid leaving. The customary index of comparison 
between actual behaviour and the theoretical, or ideal, 
behaviour is efficiency. There are several forms of plate 
efficiency, which may be related in various ways. It is 
important, therefore, to define plate efficiency, and to 
realize which of the plate efficiencies are measured under 
particular experimental conditions. 


Overall plate efficiency E~ 


One of the most useful forms of plate efficiency is due to 
Lewis,! the overall plate efficiency E°. It is defined as the 
ratio between the number of theoretical plates necessary 
to effect a given separation and the number of actual plates 
required for the same separation. This overall plate 
efficiency has, in itself, no fundamental mass transfer 
basis, but is still quite useful since only terminal conditions 
are required for its application. 


Murphree point efficiency £,., 

The Murphree point efficiency £,,, was derived by 
Murphree,'* who assumed a simplified expression for the 
instantaneous mass transfer to one bubble of vapour 
rising through the liquid. The derivation assumes that the 
vapour composition in the bubble changes continuously 
by mass transfer, and that the bubble is rising through a 
liquid of constant composition. It is defined as: 


Fe. = eee meyb 


where Y’ is the local mole fraction of the more volatile 


1 
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component in the bubble entering the liquid, Y,; is the 
local mole fraction composition of the bubble leaving the 
liquid, and Y, is the mole fraction composition of the 
vapour in equilibrium with the liquid of uniform compo- 
sition on the plate. 

Murphree" has also shown that E;,, may be expressed as 


Ee =1—exp (—m) whet) 


PK, A@ 
where m= —= 


Lake) 


and P is the total pressure, Kg is the overall mass transfer 
coefficient, A is the interfacial area of the bubble, @ is the 
contact time, and G is the moles of gas in the bubble. 

In many cases there are considerable differences in the 
composition of the liquid at various points on the plate. 
In these instances, E,,, applies as derived for any point on 
the plate, but not for the whole plate. 


Murphree plate efficiency £,,, 


The conditions necessary for the proper application of 
E‘., may not be met under the defined experimental con- 
ditions. However, it is still convenient to have a form of 
efficiency similar to E,,,, but defined for the plate as a whole. 
This form is known as the Murphree plate efficiency E,,,, 
and is defined: 


pee 2 anid) 


where Y,and Y, are the average composition of the vapour 
entering and leaving the plate, respectively, and Y= is the 
composition of the vapour in equilibrium with the liquid 
leaving the plate. 

The equations for E,,, and E,,, appear similar, but the 
former uses values for a small region of the plate, while 
the latter uses average values. In general, E,,, may be con- 
sidered the integrated effect of all the E;,, values on the 
plate. 


Relation between E° and E,, 


Lewis! has shown a relation between the overall column 
efficiency E° and the Murphree plate efficiency E,,,. This 
relation is 


Few MF Ee A= 1) F118) 
Ind 

ak, 

a 


where mn 


and K is the equilibrium constant ( Y/X), V is the vapour 
rate, and L is the liquid rate. 

This relation is correct provided that: (1) L/V is con- 
stant; (2) the equilibrium line is straight, Y, = KX+); 
and (3) E;,, is the same for all plates considered. 
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Figure 3. Overall plate efficiency E° versus Murphree plate 
efficiency Emy 


Equation (5) is presented graphically in Figure 3. It can 
be seen that as long as (KV/L) is close to 1:00, E° and E,, 
differ by only a few per cent. 

In the work with the small diameter columns, the ratio 
(KV/L) was close to 1-0, and hence there should be no 
significant difference between values of E° and E,,. The 
validity of this conclusion can be seen in Figure 4, where 
plate efficiencies of the 1-35 in. diameter columns previ- 
ously shown in Figure 2 are plotted against vapour 
velocity. The number of plates for which the efficiency was 
measured is now given as the parameter. Those values for 
one plate should correspond to E,,. Those values for two, 
three, and four plates should correspond to E°. There is 
apparently no significant difference between the several 
sets of efficiencies, indicating that E° and E,,, are approxi- 
mately the same. 


Relation between F,,, and E’.,, 


On the theoretical plate, the average composition of the 
vapour leaving is the equilibrium value for the liquid 
leaving. If the vapour and liquid ona plate were completely 
mixed, it would be impossible to obtain a better separation 
than that given by a theoretical plate, and the value of E<, 
approaches that of E,,,. However, when there is a compo- 
sition gradient in the liquid across the plate, the average 
concentration of the more volatile component on the plate 
may be greater than the concentration of the liquid 
leaving the plate. As a result of this greater concentration, 
the vapour actually leaving the plate may have a concen- 


tration greater than that in equilibrium with the liquid 
leaving the plate. 
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It is thus possible for this concentration gradient effect 
to give E° (and hence E,,) values greater than the E/, 
values, since such gradients do not apply to the E,,, values. 

Lewis?* also considered the effect of several experimental 
conditions, assuming (1) E,,, was constant across the 
plate, and (2) Y, = KX+b. For complete vapour mixing, 
but negligible liquid mixing on each plate, Lewis derives 


eres OE Leeda! ...(6) 


Figure 5 presents equation (6) graphically, with 
A = (KV /L) as the parameter. It can be seen that, under the 
condition of liquid non-mixing described, the values of 
Envy (and hence E°) can exceed the values of E”, 
significantly. 
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Figure 4. Plate efficiency; number of plates as the parameter (1-35 in. 
columns, sieve plates) 


Discussion of experimental results 


Figure 2 indicated that the plate efficiencies measured 
in the small columns were considerably larger than those 
of the pilot plant and those of Sellers and Augood. In the 
small columns, the weirs were-directly opposite one 
another, so that the liquid had a direct and unobstructed 
path. In this case, it is possible that the liquid was able to 
flow directly across the plate with little or no mixing. Hence 
the efficiencies measured would correspond to values of 
E,y (or E°) and not to £;,,. 

In the pilot plant column, the downcomers were 
adjacent to one another, separated bya weir. This arrange- 
ment caused a semicircular liquid path which tends to 
induce a more complete liquid mixing. 

In the columns of Sellers and Augood,® the liquid path 
was interrupted by a bubble cap, which undoubtedly 
caused complete liquid mixing. 

Since complete liquid mixing occurred in the columns of 
Sellers and Augood, and since their data fits that of the 
pilot plant column, the case for complete mixing in the 
pilot plant column is reinforced. In these instances, the 
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Figure 5. Plate efficiency Ef,y versus point efficiency Ejny 


efficiencies measured would correspond to values of E,,, 
Figure 6 illustrates this argument. 

In Figure 6, the best curve of the plate efficiency data 
from the small diameter column is repeated from Figure 2 
as the dashed line. The best curve of the data from the pilot 
plant column and from Sellers and Augood is repeated 
from Figure 2 as the solid line. 

The points in Figure 6 are the result of treating the data 
from the small columns by means of equations (5) and (6). 
In this treatment, it was assumed that this data actually 
represented E., (or E°), corresponding to incomplete 
liquid mixing. Then equation (6) was applied to determine 
the corresponding E,,, value for complete mixing. The data 
treated in this manner is seen to correspond quite well with 
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Augood) 
© Data of 1:35 in. column treated for liquid mixing 


Figure 6. Plate efficiencies of the 1:35 in. columns treated for liquid 
non-mixing 
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the experimental curves of the pilot plant and of Sellers 
and Augood. 

The results of this treatment, as presented in Figure 6, 
offer a strong argument for the case of incomplete mixing 
in the small columns, and for complete mixing in the pilot 
plant and in the columns of Sellers and Augood. In the 
former case, the efficiencies actually measured correspond 
to E,, (or E°), and in the latter cases to E/,. 


Table 1. Data of typical run 


Flow Rates 
Stream (ft9(n.t.p.)/min) (g.moles/min) 
Bottoms 0-11 0:14 
Feed 14-16 17:89 
Distillate 14-05 17:75 
Circulation 12:31 15:55 
Pressure, absolute 
Point (atm) 
Circulation (interior of reboiler) 4-0 
Column (exterior of reboiler) 1:61 
Hydrogen reservoir 0:83 
Temperature 
Point (°K) 
Bottoms (exterior of reboiler) 24:5 
Feed (into the column) 39-7 
Distillate (top of column) 22:1 
Circulation 
(into reboiler) 28-0 
(interior of reboiler) 26:0 
Hydrogen reservoir 19-8 
Analyses 
Composition 
Sample (mole % D,) 
Bottoms 96:9 
Feed 0:77 
Distillate 0-015 
Circulation 0-015 
Between plates 28/29 87:0 
Between plates 13/14 1-:0— 
Between plates 11/12 1-:0— 
Between plates 2/3 0-035 
Conclusions 


A few conclusions may be inferred from this analysis. 
The first is that plate efficiency is not per se a function of 
the column diameter, as might have been supposed. 
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Instead, the plate efficiency depends upon the degree of 
liquid mixing. If this mixing is complete, the efficiency is 
invariant with respect to plate diameter. 

On the small diameter perforated plates, it is possible 
to have incomplete mixing, corresponding to a higher 
value of apparent efficiency. 

It should be possible to obtain higher values of plate 
efficiencies in hydrogen distillation columns by inhibiting 
the liquid mixing, since the E° (or E,,,) values are higher 
than the E7., values. The feasibility of this suggestion 
depends at least upon the arrangement of the weirs and 
downcomers, the column diameter, the vapour and liquid 
rates, and the mass transfer liquid mixing length. Never- 
theless, this suggestion may prove attractive. 

The comparison of efficiencies indicates that efficiency 
does not depend strongly upon the type of plate. Since the 
efficiency of perforated and bubble cap plates is much the 
same, the choice of plate depends more strongly on the 
plate hydraulic characteristics. 
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IN 1932 Keesom and Kok! investigated the specific heat 
anomaly at the superconductive transition point of tin. 
Two of their values, however, differed widely from the rest 
and this was ascribed to the accidental presence of a small 
magnetic field. On the basis of this observation, Mendels- 
sohn suggested that the adiabatic transition between the 
superconductive and the normal state might be accom- 
panied by a magneto-caloric effect which could be used for 
obtaining low temperatures. Magnetization experiments 
on tin by Mendelssohn and Moore? demonstrated the 
existence of this effect and small coolings (~ 0-3° K) were, 
in fact, observed. This result was confirmed by Keesom 
and Kok? who obtained coolings of the order of 0-15°K in 
thallium. A somewhat improved arrangement* yielded 
final temperatures of the order of 0-2°K which were 
estimated from an extrapolation of the critical field curve 
of tin to low temperatures. It was decided, prior to the 
development of the cooling method, to measure more 
carefully the critical fields of a number of superconductors 
which might serve as cooling agents. This work® was inter- 
rupted by the war. Work on adiabatic magnetization was 
not resumed until 1954 when the present series of experi- 
ments was started. 


Thermodynamic considerations 


Allsuperconductors below their transition temperatures 
can be brought back to their normal state by the applica- 
tion of a strong enough magnetic field called the critical 
field H,. H, rises from zero at the transition temperature 
to a temperature independent value as we approach the 
absolute zero. Since the thermodynamics of the transition 
is analogous to any other phase transition, one can calcu- 
late the entropy differences in the two phases as 


ss = ed He 
Stet ‘ot iegenid Tl bo" UP 


where S, and S, are the entropies in the normal and 
superconducting states, respectively. Since dH,/dT is 
always negative, we have 


ue Ge 


except when H, = 0, or dH,,/dT = 0, ie. at the transition 
temperature 7, and at the absolute zero where dH,/dT 
must vanish as a consequence of the third law of thermo- 
dynamics. 
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The superconducting state, therefore, is a more ordered 
state of the metal and is due to a fraction of the conduction 
electrons going into a state of zero entropy below the 
transition temperature.® The temperature—entropy curve, 
which is typical for a superconductor in both phases, is 
given in Figure 1. There are two electronic contributions 
S, and S, corresponding to the normal and supercon- 
ductive states, and the contribution S),, due to the lattice 
vibrations, which is common to both. 


Figure 1. Diagram of entropy S versus temperature T curve 


From equation (1) it is clear that in isothermal des- 
truction of superconductivity by a magnetic field, heat 
will be absorbed, and, if the magnetization is carried out 
adiabatically, the temperature must fall. Thus a super- 
conductor cooled in zero field to P will, when magnetized 


. adiabatically, cool until it reaches Q where its entropy in 


the normal state is the same as in the superconducting 
state at P. 

Corak, Goodman, Satterthwaite, and Wexler’ and 
Goodman? have shown from calorimetric measurements 
that, well below the transition temperature, the electronic 
contribution to the atomic heat in the superconducting 
state can be represented by 


yy 
Cry Syriaep beg ETD) 


where y is the electronic specific heat coefficient in the 
normal state, T, the transition temperature, and a and 5 
are constants depending on the metal. The above expres- 
sion has also been derived theoretically by Bardeen, 
Cooper, and Schrieffer? suggesting a ~ 8-5 and b ~ 1:44 
for all superconducting metals. From equation (2) we can 
calculate the cooling ideally obtainable by equating the 
entropies in the two phases. The relation between the 
initial temperature 7, and the final temperature 7; obtained 
on magnetization is 
WR ee 3 3 
Tea | exp(~6T./T) 47, A T; A T; 
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where () is the Debye characteristic temperature of the 
metal and A has the value 464-4. In this way, by calculating 
7, for various values of T; we have plotted curves showing 
the relation between the initial and final temperatures for 
tantalum and tin (Figure 2). From these curves we see that 
starting from 1-25°K, final temperatures of 0-05°K and 
0-1°K can be attained on magnetizing tantalum and tin, 
respectively. 


(°K) 


|;— 


—— 
0) 0:5 1-0 1-5 2:0 ZS 
7; — (°K) 


Figure 2. Diagram showing the relation between the initial and final 
temperatures for tantalum and tin 


Advantages and disadvantages of the method 


Comparing this method with the adiabatic demagnet- 
ization of paramagnetic salts, we find that it has certain 
disadvantages. The most important of these is the linear 
fall with temperature of the specific heat leading to low 
absolute values near 0-1°K. This makes the cooling 
efficiency in this region, measured in calories per cubic 
centimetre, much smaller than that for paramagnetic 
salts, and although fairly large temperature drops may be 
effected, it will be difficult to maintain the lowest obtain- 
able temperatures for any length of time. The smallness of 
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the entropy changes near 0-1°K also makes this process 
relatively inefficient for cooling appreciable quantities of 
other substances. 

Another aspect which tends to reduce the efficiency of 
this method is the production of irreversible Joule heat. 
During the process of magnetization the specimen 
traverses the intermediate state and a change of field with 
time induces eddy currents in the normal regions whose 
dying out causes heating of the specimen. Consequently, 
unless the rate of increase of the field is very slow, the 
temperature drop may be considerably less than that 
calculated from equation (3). 

There are, however, several points in favour of this 
method. For instance, the necessary fields to destroy super- 
conductivity are of the order of a kilo-oersted or less for 
all superconductors and can be easily produced over large 
volumes. One can, therefore, magnetize considerable 
quantities of a metal and in this way compensate, to some 
extent, for the smallness of the heat capacity. Another 
attractive feature of the method is the comparatively large 
thermal conductivity of metals below 1°K.?° Near 0:05°K 
metals have a conductivity twenty times larger than that of 
paramagnetic salts.1! This fact coupled with the small heat 
capacity makes the equilibrium times very small. The 
presence of a metal also makes it easier to make thermal 
contacts, and specimens to be investigated can be directly 
soldered to the refrigerant metal. 

It appears that the method would be particularly useful 
in the temperature range between 0-3 and 1°K. Most 
paramagnetic salts in common use have specific heat 
maxima below 0-3°K and their heat capacities between 
0-5 and 1°K are rather small. On account of the fact that, 
due to gas desorption, vacuum conditions begin to de- 
teriorate above 0-6°K, work with them becomes increas- 
ingly difficult above this temperature. Moreover, the 
specific heat of the salts fall with rise in temperature and 
this accelerates the warming-up rates. In the supercon- 
ducting method, on the other hand, the entropy of the 
final phase is rising linearly with absolute temperature and 
the rate of warming-up gets retarded by the steadily 
increasing heat capacity. 


The choice of the metal 


Metals. with a high Debye ©, a high transition tem- 
perature, and a large coefficient of electronic specific heat 
y, are particularly suitable as cooling agents. Since the 
amount of heat which has to be extracted from the metal 
to cool it from hydrogen to helium temperatures is 
inversely proportional to the third power of © and, since 
a fairly large quantity of the metal has to be used, a high 
value of © is desirable in order to reduce the quantity of 
liquid helium required for preliminary cooling. Tantalum 
and niobium satisfy these conditions, but the quantities 
required to cool a reasonable amount of any other sub- 
stance (say 10 cm? of liquid helium) to an appreciable 
degree is of the order of 25 moles, which was considered 
excessive for our purposes. Tin, on the other hand, is 
readily available in large quantities at a reasonable cost. 


CRYOGENICS - DECEMBER 1960 


It was therefore decided to use tin, although it has a much 
smaller Debye © and, on account also of a smaller y than 
that of tantalum and niobium, is not as good a cooling 
agent. 


Apparatus 


The apparatus had to be designed with the following 
essential features: 


(1) A large cryostat capable of accommodating and cool- 
ing about 40 moles of tin to a temperature of about 
1K; a 

(2) A helium vessel surrounding this cryostat maintained 
at a temperature of about 1°K, so as to reduce the 
heat influx; and 

(3) A suitable thermometer to measure accurately the 
temperature in the range of about 0-1 to 1°K. 


Condensing line 
for helium 


Thin stainless 
steel wire 


Vacuum 


Figure 3. The apparatus 


The apparatus in its final form is shown diagram- 
matically in Figure 3. It consists of a cylindrical brass 
cryostat 30cm long and 8cm in diameter which can hold a 
block of tin weighing 4-5 kg. The surrounding helium 
vessel with a total capacity of 800 cm® was made up of two 
separate parts. The upper half consists of a copper vessel 
of 700 cm? and the lower is made up of two concentric 
brass cylinders, the space between which contains the 
remaining 100 cm? of helium. A composite helium vessel 
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was necessary in order to facilitate the suspension and 
removal of the tin block which is held in position inside the 
cryostat by means of thin stainless steel wires. A small 
vessel V, made of thin walled copper, having a capacity of 
about 2 cm? is soldered on the top of the tin specimen. By 
condensing a small quantity of liquid helium into this 
vessel and pumping it off vigorously, the temperature of 
the specimen could be lowered to 0:8°K before magnetiz- 
ation, while the outer helium bath was at 1-:25°K. This 
procedure was adopted after the failure of a number of 
experiments, in which the specimen was cooled to 1:25°K 
by means of helium exchange gas. Then, magnetization 
after the vacuum space had been pumped for several hours 
resulted in no cooling at all, showing that the exchange gas 
could not be successfully pumped out at 1-25°K. More- 
over, by pumping the helium in the vessel V the specimen 
could be magnetized from a lower initial temperature. 


Temperature measurement 


For measuring the temperature of the specimen we used 
0-25 W LAB resistors T, and T, (Figure 3) with a nominal 
value of 33 Q at room temperature. Their resistance varied 
from about 160 Q at 4:2°K to about 2,000 2 at 1-0°K. 
while their sensitivity defined as dR/RdT rose from 
approximately 0-5 at 1:0°K to above 10 at 0:2°K. A 24 
s.w.g. enamelled copper wire was tightly wound round 
the carbon and a heavy layer of ‘Araldite’ was applied all 
round it to cover completely the body of the resistor. It 
was then baked until the ‘Araldite’ was set hard, giving a 
good thermal contact between the copper wire and the 
carbon. The loose ends of the copper wire were then 
soldered to the tin specimen. Thermometer resistance was 
measured by a Wheatstone bridge having equal ratio 
arms and, to cancel the lead resistance, a three-lead con- 
nection to the resistor was used. 


Calibration of carbon resistance thermometers 


The lowest temperature to which we could calibrate our 
thermometers with the above apparatus was 1:25°K 
measured by means of helium vapour pressure mano- 
meters. However, the range it was expected to cover by 
magnetizing tin was between 0-1 and 1°K and, therefore, 
calibration against a susceptibility thermometer below 
1:25°K was necessary. In the early experiments the re- 
sistors were calibrated in a separate demagnetization 
cryostat and then attached to the specimen. It was found 
that the process of transferring from one cryostat to the 
other, which involved soldering, invariably changed their 
characteristics by about 3 per cent and in some cases errors 
of up to 15 per cent were introduced in temperature de- 
terminations. If, however, the resistors were left untouched 
in between experiments, their reproducibility remained 
within | or 2 per cent. Since the same arrangement was 
used for specific heat measurements, several experiments 
were made to see under what conditions maximum repro- 
ducibility could be achieved. It was found that if the 
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following precautions were observed the thermometers 
were reproducible to better than 0:1 per cent: 


(1) The carbon of the thermometer was completely 
covered in ‘Araldite’ before baking, thus preventing 
moisture from condensing on the carbon during the 
experiments ; 

(2) It was soldered to the specimen once and for all and 
left undisturbed from one experiment to another; and 

(3) The thermometers had to be cooled to helium tem- 
peratures from room temperature three or four times 
before calibrating. 


Hence, it was necessary to calibrate the resistors against 
magnetic susceptibility in contact with the tin specimen 
so that they were left untouched between a calibration and 
a magnetization experiment. For this purpose an alterna- 
tive set of vacuum cases and a helium vessel were designed 
which had appropriate tails at the bottom to accommodate 
the salt pill and the measuring coil. By demagnetizing the 
salt pill the temperature of the specimen could be lowered 
to 0-15°K and the thermometer resistance calibrated 
against the susceptibility of the salt. The cryostat with 
the salt pill in position and in thermal contact with the tin 
specimen is shown diagrammatically in Figure 4. 
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Figure 4. The cryostat with the salt pill in position 
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The magnetic field 


For magnetizing the specimen a transverse field pro- 
duced by a horseshoe type of electromagnet was used. It 
was capable of producing a field of about 1 kgauss. In the 
earlier experiments it was excited by a battery source and 
in the later, for technical reasons, by d.c. mains. 


Experimental results and discussion 


The specimen used in the early experiments was a 
cylindrical block of tin 15cm long and 6-95 cm in diameter. 
It was supplied by the Tin Research Institute, London, and 
its purity as quoted by the suppliers was 99-990 + 0-002 
per cent. It was annealed at a temperature of 200°C for 
2 weeks before magnetization. 

Theoretically, the transition can be made isentropic 
only if the field is applied infinitely slowly. For practical 
purposes an order of magnitude calculation showed that, 
in order to keep the gain in entropy due to eddy current 
heating negligible, the magnetization would have to take 
several minutes. Since owing to the complicated geometry 
no exact calculations were possible, it was decided to study 
the effect of varying the rate of increase of the field. 
Results of the coolings obtained on applying a field of 305 
gauss at various rates are given in Table 1. 


Table 1 


Time taken to apply _ Final temperature 
the field achieved 
(min) (°K) 


2 0-6 
4 0-52 
6 0-44 
8 0-42 

10 0-38 

12 0-36 

14 0:35 

16 0:35 


Increasing the time of application of the field beyond 
16 min did not result in reaching a lower temperature, and 
0-35°K was, in fact, the lowest obtained in this series of 
experiments. This is a good deal higher than the expected 
value of 0-02°K and shows that the magnetic heating 
caused by inducing the transition in a finite time is not the 
only factor responsible for reducing the temperature 
changes achieved on magnetization. Dolecek in his at- 
tempts to study the temperature—entropy relations of lead 
spheres has also found that the fall in temperature on 
adiabatic magnetization was considerably smaller, pro- 
vided the temperature before magnetization was well 
below T.. The reduced change in temperature could be 
caused by the following reasons. 

(1) If, as is possible in a polycrystalline specimen of large 
dimensions, there exist regions of strain with high critical 
fields, then a fraction of the metal may remain supercon- 
ducting and cause a higher final temperature. But the 
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fraction which would give a final temperature of 0-35°K 
is nearly 0-8 and is an impossibly high value. However, in 
order to make sure, experiments were repeated after 
annealing the specimen for a further period of 2 weeks at 
a temperature of 200°C, but no improvement was found. 

(2) A wellannealed specimen of pure tin is likely to have 
large flawless crystals which may superheat and cause eddy 
current heating independent of the rate of change of the 
field. Although no accurate measurements were made, 
evidence of superheating at the beginning of the transition 
was found on magnetizing an ellipsoid of pure tin with a 
known demagnetizing coefficient. Just before the field 
began to penetrate, the value of the field at the equator of 
the ellipsoid was 3 to 4 oersteds greater than the critical 
field H, for the initial temperature. 

(3) The process of increasing the field by reducing the 
resistance in the circuit of the electromagnet could not be 
made continuous by using an ordinary rheostat which 
changes its resistance in steps. Due to this, there would be 
a certain amount of Joule heating independent of the time 
taken to complete the transition. In order to reduce the 
heating due to this cause, subsequent experiments were 
performed on a laminated tin specimen which consisted 
of 300 circular disks of pure tin 0-5 mm thick and 6 cm in 
diameter. They were tightly packed in a brass cylinder and 
were insulated from each other by means of a thin layer of 
varnish. In order to improve thermal contact between the 
disks and the cylinder, the dead space was filled with 
helium at a pressure of 15 atm at room temperature. 
Magnetization of this composite specimen gave final 
temperatures considerably higher than 0-35°K. Besides, 
the time taken for preliminary cooling to 4-2° K was over 
50 per cent longer than in the case of a solid cylinder. It 
was also noticed that the performance of this specimen 
deteriorated with each subsequent experiment. This was 
later found to be due to a gradual transformation of white 
tin to the grey form. On account of these difficulties, 
experiments with the laminated specimen were abandoned 
and it was decided to reduce the eddy current heating by 
alloying the specimen; this increases the resistivity of the 
metal and also prevents superheating. In choosing the 
impurity for tin, the following considerations were 
important: 


(1) The impurity should form a solution in tin so that the 
alloy is as strain-free as possible. This was necessary 
to prevent local regions of very high critical fields due 
to the so-called ‘Mendelssohn sponge’ ;}” and 

(2) The thermodynamic properties of tin should not be 
altered drastically by the impurity. 


Indium, which has the same order of atomic diameter as 
tin, is known to form a substitutional solution with tin in 
concentrations up to 3 per cent. Doidge!* has shown that 
the difference in thermodynamic properties of tin and 
tin—indium alloys up to 3 per cent is very small. Transition 
temperatures are only slightly lower than that of tin and 
the critical field curves indicate that the electronic specific 
heat remains virtually unchanged. 
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In order to study the effect of alloying, three specimens 
containing 0-5, 1, and 2 per cent indium were prepared. 
Their sizes and shapes were the same as the original 
cylinder of pure tin. The use of alloys with a higher con- 
centration of indium was not contemplated because of the 
evidence?* that the interphase energy becomes negative at 
about 2:5 per cent of indium. According to H. London,"4 a 
negative interphase energy can cause thin superconducting 
filaments to persist in fields greater than the critical, and it 
was feared that these might cause a reduced drop in tem- 
perature on magnetization. All the three specimens were 
annealed at a temperature of 200°C for 3 weeks before 
magnetization. The lowest temperatures obtained are 
shown in Table 2. 


Table 2 


Indium concentration Final temperature 


(per cent (wt)) CK) 
0:5 0-192 
1:0 0-180 
2:0 0:345 


The temperatures 0-192° K and 0-180°K, although sub- 
stantially lower than 0-35°K attained with pure tin, are 
still a good deal higher than the expected value of 0-02° K. 
These specimens were magnetized by using d.c. mains as 
a source of excitation for the electromagnet, and this may 
have possibly resulted in an excessive heating during the 
magnetization. The final temperature of 0-345° K obtained 
with the 2 per cent specimen was due to the Mendelssohn 
sponge behaviour. This was evident from the fact that on 
increasing the field beyond 305 oersteds the specimen 
showed slight further cooling in small steps. This, however, 
could not be utilized for obtaining a lower temperature 
because in raising the field the specimen warmed up as a 
whole. 
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Figure 5. Diagram of rise in temperature versus time curve for pure 
tin in a steady field 
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Warming-up rates 


In order to assess quantitatively the value of the method, 
it Was necessary to measure the rate at which the specimens 
warm up after magnetization. Figure 5, which is a plot of 
the rise in temperature against time for pure tin in a steady 
(battery) field of 305 oersteds, shows that it takes more than 
6 hr to reach 0-75°K. This compares favourably with the 
warming-up rates of an average demagnetization cryostat 
in the same range. The alloys which were magnetized in 
the fluctuating mains field showed a much faster rate of 
warming, as can be seen from Figure 6. On a more careful 
examination of the curves in Figure 6, we find that the 
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Figure 6. Diagram of temperature versus time curves for pure tin and 
tin-indium alloys : 


1 per cent alloy has the fastest average rate of rise in 
temperature and the 0:5 per cent the slowest. This seems 
surprising because, if the increase in the rate of warming- 
up is caused by Joule heating in the normal metal the 
2 per cent alloy should have the slowest rate and the 0-5 
per cent the fastest, since Joule heating is inversely pro- 
portional to the resistivity of the metal. The frequency of 
the ripple and the resistivities of our alloys were such that 
skin depth considerations were unimportant. Another 
feature of these curves is that the 0-5 percentalloy, although 
it warms up more slowly on the average, between 0-2 
and 0-4°K its rate of rise in temperature is greater than 
that for the 1 per cent alloy, which warms up at a higher 
rate on the whole. At 0-4°K, a temperature which is 
common to all the three specimens, the slopes of these 
curves are approximately in the following ratios: 


O55 ue eo ee 1 t-37.0-8 
whereas at 0-6° K they change over to 
0:5% :1%:2% = 1-4:0-8:1-0 


This behaviour could be explained either by assuming 
that the resistivities of these alloys were changing with 
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temperature, or by a change of specific heats, or by both. 
The likelihood of a change in specific heat was rather small, 
since a fairly rapid change would have to be assumed to 
explain the result. On the other hand, a change of effective 
resistance with temperature in the presence of a magnetic 
field is possible. The fields required to restore completely 
the resistance in an alloy are much higher than the fields 
necessary for the flux to penetrate. This means that a very 
small fraction of the total volume of the alloy remains 
superconducting and has a higher critical field curve. Thus, 
it is possible that because of field fluctuations a portion of 
the material at a particular temperature alternates between 
the normal and the superconducting phase, causing a good 
deal of irreversible heating. The rate of warming-up at a 
particular temperature in alloys would depend on the 
amount of material which alternates between the two 
phases and on the residual resistance of this material. On 
warming-up, more and more of this material becomes 
permanently normal, until a temperature is reached where 
the transition to the normal state is complete and the 
warming-up becomes steady. Beyond this point the heating 
rate will depend only on the resistivity of the alloy. 

This explanation is supported by a more detailed 
investigation of the warming-up curves over a wider 
temperature range. In the 0-5 and 1 per cent alloys it was 
noticed that the rate of warming-up showed discontinuous 
changes. Assuming that this was due to a sudden change in 
heat input and not to a discontinuous change in specific 
heat, the curves shown in Figure 7 were plotted. Here, the 
specific heat for all three specimens was supposed to be 
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Figure 7. Diagram of energy input versus temperature curves for pure 
tin and tin-indium alloys 


CRYOGENICS - DECEMBER 1966 


the same as for pure tin. Using the empirical values of 
Corak and Satterthwaite,! the energy input per minute 
was calculated from the observations and plotted against 
temperature. All the three curves show marked discon- 
tinuities, which can be explained by assuming that at 
certain temperatures a small but finite amount of material 
goes over to the normal state permanently and causes a 
corresponding decrease in irreversible heating. This goes 
on until the whole of the alloy is normal. The curves then 
cease to show these discontinuities and attain a steady 
value. For the 0:5 per cent alloy this seems to occur at 
0-7°K and for the | per cent alloy at about 0-88°K. The 
2 per cent alloy, which is likely to have greater inhomo- 
geneities, does not seem to lose all its superconducting 
material until about 1:5° K. Above this temperature, where 
all the specimens show a steady heat input, we find that 
the lowest value is for the 2 per cent alloy and the highest 
for the 0-5 per cent alloy. This seems eminently reasonable, 
because their electrical conductivities are also in the same 
sequence. 

In spite of the excessive heating caused by the ripple 
field, the rate of temperature rise in our specimens was 
slow enough to allow quantitative measurements below 
1°K. This was proved in subsequent experiments by mak- 
ing measurements of specific heats of our specimens down 
to 0:-5°K. Moreover, the mechanism suggested here is 
further supported by our measurements of the specific 
heat (to be reported elsewhere). 


Low temperature thermometry and its accuracy 


The errors involved in measuring the low temperatures 
achieved depended mainly on the sensitivity of the mag- 
netic thermometer used for calibrating the resistors. It 
consisted of a conventional type of mutual inductance 
bridge using an external compensating coil. On switching 
offa current of 100 mA in the primary circuit of this bridge, 
a galvanometer deflection of 365 mm was obtained at 1°K 
on a scale which was 4 m distant. The galvanometer de- 
flection for all measurements below 1° K was kept roughly 
at 400 mm and a tané correction of 0-5 per cent was 
applied. The slope of the deflection against 1/(T7— A) line, 
where A is the Curie-Weiss constant, was calculated by 
the method of least squares and gave a value of 364:4+0-9, 
which was used for calculating the temperature below the 
helium range. Since the galvanometer deflection could be 
measured to about one part in 1,500, the errors in the 
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temperatures were governed by the probable error in the 
slope quoted above. A larger probable error than the 
accuracy of measurement was caused by the fact that 
above 1°K the actual deflections were much less than 400 
mm and, consequently, the accuracy was not so great. 
Resistance measurements could be made to at least an 
accuracy of one in 4,000 and, therefore, the actual tem- 
peratures quoted are again accurate to the order obtained 
by the magnetic thermometer. We must not forget, how- 
ever, that below 0-5°K the chief error in the values of 
absolute temperature is due to the errors in the 7-7* 
relation which are at least of the order of 2 per cent. 


The author wishes to express his thanks to Dr. K. Men- 
delssohn, F.R.S., for suggesting the problem and for his 
help and encouragement. He is also grateful to the Tin 
Research Institute, London, for supplying the vast quanti- 
ties of pure tin used in the experiments. A part of this 
work was done during the tenure of Twisden Studentship 
in Experimental Physics, awarded by Trinity College, 
Cambridge. 
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EXPERIMENTS on unsaturated helium films have shown that 
phenomena associated with superfluidity in liquid helium 
are depressed to lower temperatures when the thickness 1s 
reduced below that of the saturated film (~ 10-§ cm). For 
instance, the work of Bowers, Brewer, and Mendelssohn?” 
and of Long and Meyer *:4 demonstrated that the tempera- 
ture of onset of superflow is lower in such thin films than 
the normal lambda temperature. Frederikse** and Aston 
and Mastrangelo’ showed that the temperature of the 
specific heat maximum also decreases with film thickness. 
Comparison of these two types of experiment reveals an 
unexpected feature, namely that the specific heat maxima 
are consistently higher than the onset temperatures for 
superflow by about 0-4 K, while in the bulk liquid both 
of the temperatures are the same within the experimental 
accuracy. 

A natural extension of these experiments is the investi- 
gation of ‘bulk’ liquid helium confined in channels of 
dimensions similar to those of the unsaturated film. 
‘Vycor’ porous glass. described by Nordberg® provides a 
suitable material for such experiments, having a sponge- 
like structure permeated by pores whose mean diameter 
depends upon the sample used, and varies between about 
30 and 100 A. Atkins. Condon, and Seki? showed that 
the temperature of onset of superflow through specimens 
of this glass was indeed depressed below the ordinary 
lambda point. 

In some respects. experiments using ‘ Vycor’ are more 
difficult to interpret than those on unsaturated films, since 
the pore size distribution is not known accurately, and the 
analysis of results is thus complicated. On the other hand, 
when comparing different aspects of the lambda transition 
with one another, ‘Vycor’ provides the advantage that 
for both types of experiment the pore size is the same, so 
that flow and specific heat measurements can be compared 
with a good degree of reliability, whereas with unsaturated 
films, the phenomena must be correlated with film thick- 
ness through the adsorption isotherms, which have been 
found to give inconsistent results.!° 1 With this situation in 
mind, we have measured the flow properties and specific 
heat of liquid helium in ‘Vycor’ glass as a function of 
temperature using the same sample for both sets of 
experiments. The pore diameter has been measured in the 
course of a separate investigation of adsorption, again in 
the same specimens, by Brewer and Champeney.!2 A brief 
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report of some of the results in the present paper has 
appeared elsewhere."* 


Flow experiments 


Experiments were carried out to determine the rate of 
flow at a given pressure head as a function of temperature 
for several specimens of ‘ Vycor’ glass. From these meas- 
urements, the temperatures of onset of superflow were 
determined, and suitable specimens chosen for use in 
specific heat measurements. The principle of the method 
was that generally used for isothermal gravitational flow 
experiments on liquid helium. A glass reservoir was first 
filled by total immersion in a bath of liquid helium, then 
partially lifted out and the level within it measured as a 
function of time while liquid flowed out through a disk of 
‘Vycor’ sealed to the bottom of the vessel. The flow rate 
at a given pressure head was determined from the graph, 
and the experiment repeated at different temperatures. 

A typical form of the reservoir is shown in Figure 1. 
The side tube was a capillary of 1 mm bore which provided 
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Figure 1. The flow apparatus showing ‘ Vycor’ disk and clamp 
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a clear meniscus for viewing. A similar tube was fixed 
beside it, open at both ends and dipping into the bath, for 
observing the bath level and compensating for capillary 
rise in the side tube. A scale engraved in half millimetres 
was supported just behind these tubes, and the levels esti- 
mated to a tenth of a millimetre by observation through 
the glass Dewar system using a long focus microscope. 

The ‘Vycor’ disk, of thickness 5 mm and. diameter 
3-6 cm, was sealed to the reservoir using a system involving 
indium gaskets. The disk (A) was sandwiched between 
two brass plates (B), with indium gaskets (C) 1-5 mm 
thick above and beloW, fitting into grooves in the brass. 
The assembly was tightened by means of six bolts around 
the rim with spring washers (D), and was connected to the 
glass reservoir with a copper-to-glass seal. 

The area of the disk accessible to flow was limited by the 
diameter of the upper gasket, the lower one acting solely 
as a cushion to prevent the disk from cracking. Two sizes 
of gasket could be used having nominal diameters, before 
compression, of 2 cm and 0-5 cm. Under pressure they 
flattened out, and the actual diameter effective in describ- 
ing the area for liquid flow was different from the nominal 
diameter. 

Experiments were performed with three reservoirs of 
different dimensions so that a large range of flow rates 
could be observed without the rate of fall of meniscus 
becoming too quick to measure. The gasket diameter d,, 
reservoir diameter d,, and reservoir orifice diameter d, for 
the three flow vessels are given in Table 1. The gasket 
diameters refer to the measured internal diameters after 
compression. The orifice diameter is included because of 
its importance in the effective action of the distillation 
process discussed below, which maintains the reservoir 


and bath at similar temperatures. Flow vessels I, II, and 
IIT were used, respectively, in the temperature ranges 1-22 
to 1-97°K, 1-84 to 2:05°K, and 2:04 to 2:11°K, and the 
results obtained were pieced together to give information 
over the whole range. 

Other ‘Vycor’ specimens were available in the form of 
‘thimbles’ or short test-tubes. For these, a simpler form 
of apparatus was used, consisting of a glass reservoir of 
suitable diameter, constricted at the top to limit film flow, 
and attached to the ‘ Vycor’ by a graded seal. 


Table | 
dg dr do 
Flow vessel (em) (cm) (cm) 
0:35 25 0:5 
il 1-74 2°5 ue 
ve 1-74 0-35 0-11 


In Figure 2 we show typical graphs of the rate of change 
of reservoir and bath levels as a function of time at three 
temperatures for one of the ‘ Vycor’ disks. Similar curves 
were found for other specimens, and from such graphs as 
these the net rate of outflow at a given pressure head was 
determined by drawing tangents to the ‘reservoir’ curves. 
In the case of flow vessel I the flow rate was then expressed 
as cubic centimetres per minute per square centimetre of 
‘Vycor’, using the known values of d, and d,. Results 


» Reservoir 


(cm) 


(sec) ——> 


(a) 


(min) ——> 


(c) 


Figure 2. Variation of liquid levels in reservoir and bath with time for different reservoirs and temperatures; (a) T=1-22° K, reservoir I; 
(6) T=1:911° K, reservoir IT; (c) T=2:05° K, reservoir III 
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Figure 3(a). Flow rate F versus temperature for constant level 
difference of 2.cm 


similarly calculated using flow vessel II consistently 
differed from those obtained with vessel I at the same 
temperature, being four times smaller. It is reasonable to 
assume that this is due to the differing geometries of fiow 
path as in neither case would the flow have been uniformly 
perpendicular to the plane of the disk. Results with flow 
vessels II and III were therefore scaled by the same factor 
to join on to those with vessel I. This uncertainty in the 
absolute flow rate per unit area of ‘Vycor’ should not, 
however, affect its general temperature variation. 

Appreciable corrections had to be made to the flow 
rates thus determined, due to film flow and to evaporation 
from the reservoir. The effect of film flow was estimated 
from the orifice diameter and the values of film transport 
over clean glass surfaces given by Daunt and Mendels- 
sohn.'* In vessel III this correction amounted to 20 per 
cent at 2:11°K, falling to 1 per cent at 2-05°K. Otherwise 
it was negligible. 

Excess evaporation from the reservoir occurs as a result 
of the mechano-caloric effect!® when superfluid flows out 
through the ‘Vycor’ disk, and through the film. It was 
allowed for by dividing the observed rate of fall of level by 
(1+ ST/L) where Sand Lare the specific entropy and latent 
heat of helium at the temperature 7. This expression is 
based on the assumptions that no entropy is carried out of 
the reservoir by the liquid helium flowing through the 
“Vycor’ and the film, and that the effective mechano- 
caloric heat input is entirely removed by the distillation 
process discussed in detail by Atkins.!® The first of these 
assumptions is contained in London’s equation?’ for the 
thermomechanical effect, which has been tested frequently 
and verified for even larger channels than those involved 
here. The effectiveness of the distillation process was 
tested separately for each reservoir by performing repeated 
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Figure 3(b). log F versus temperature for constant level difference 
of 2.cm 


flow experiments at the same bath temperature, but with 
various currents applied to an electric heater inside the 
reservoir. Heat inputs of up to 1-5 times the effective 
mechano-caloric heat input were applied and were found 
to result simply in an increased rate of fall of meniscus 
equal to that calculated as resulting from an increased rate 
of evaporation. This fact was taken as satisfactory evidence 
that the flow rate through the ‘ Vycor’ was not affected by 
any small temperature differences between reservoir and 
bath, and also that any heat produced in the reservoir was 
removed by distillation rather than by conduction through 
the glass. Tests of this nature were performed throughout 
the temperature ranges used with each vessel. 

In order to plot the flow rate through a given ‘ Vycor’ 
specimen as a function of temperature, it is necessary to 
choose a standard pressure head for all temperatures, 
since, as the graphs of Figure 2 show, the flow is strongly 
pressure dependent. Choosing a pressure head of 2 cm, we 
obtain curves of flow rate F versus temperature as shown in 
Figure 3(a) for two specimens, in which the corrections 
described above have been included. Curves of similar 
shape are obtained for other pressure heads, but with 
different magnitudes of the flow rate. Onset temperatures 
for superflow may conveniently be obtained from the 
logarithmic plots of Figure 3(b) by a method explained in 
the final section. Qualitatively, it may be seen from 
Figure 3(a) that onset temperatures between about 1:4°K 
and 2:1°K are found in the specimens investigated. The 
lower onset temperatures were, however, obtained with 
the ‘thimble’ specimens which were not suitable for 
specific heat and adsorption measurements. For the latter 
experiments, the two disks giving curves labelled X and Y 
in Figure 3(b) were finally chosen, having onset tempera- 
tures very close to each other. 
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Specific heat experiments 


The heat capacity of helium contained in the pores of 
the ‘Vycor’ disks under saturation conditions was meas- 
ured by the usual method, employing an adiabatic vacuum 
calorimeter. The results were calculated from the tempera- 
ture increases produced by known electrical heat inputs. 
The apparatus is shown in Figure 4. 

The disks fitted, with little dead space, into a copper 
calorimeter which was suspended inside a vacuum jacket. 
The outside of the jacket was maintained at 4:2°K. The 
calorimeter had attached to it an auxiliary container into 
which a small quantity of liquid helium could be con- 
densed. By pumping on this helium, the calorimeter could 
be cooled to 1-2°K; measurements were begun after this 
liquid had completely evaporated. 

The calorimeter was filled with helium gas from a stor- 
age vessel of known volume and the amount of liquid 
helium in the calorimeter was calculated from the decrease 
in pressure of this vessel. Corrections were applied for the 
amount of gas filling the dead space in the inlet tubes. The 
same apparatus was used for measuring adsorption iso- 
therms of helium and nitrogen on these specimens of 
“Vycor’, and the volumetric system and method of allow- 
ing for dead space corrections will be described more fully 
along with that work.!” 

To reduce heat inputs the inlet tubes were coiled in 
passing from the vacuum jacket to calorimeter, and con- 
sisted of 30 cm of 1 mm diameter german silver tubing. 
The heat leak into the calorimeter was chiefly due to film 
flow along these tubes and amounted to 1-4 10-4 W at 
1-4°K, falling to a steady value of 8 x 10-§ W above 2-18°K. 

The temperature of the calorimeter was measured with 
a carbon resistance thermometer calibrated against the 
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Figure 4. Specific heat apparatus 
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Figure 5. Change in resistance thermometer versus time during a 
typical warm-up period 


vapour pressure of helium within the calorimeter, using 
the T;, experimental temperature scale and correcting for 
thermo-molecular effects at the lower temperatures as 
described by Keesom.}8 Changes in conductance of the 
thermometer with temperature, d(1/R)/dT, were measured 
to an estimated accuracy of | percent. The heater consisted 
of approximately 450 Q of 40 s.w.g. constantan wire 
wound on to the body of the calorimeter, arranged with 
potentiometer leads so that its resistance could be accur- 
ately determined. 

The heat capacity C of the system was calculated from 
the relation 

rt d(1/R) 
AACR AF 

where A(1/R) was the change in conductance of the ther- 
mometer resulting from the passage of a current / for time 
t through the heater of resistance r. In determining the 
true specific heat of the liquid helium within the *‘ Vycor’ 
pores, the following factors and corrections have been 
considered: (1) the effect of non-equilibrium within the 
calorimeter and the variation of warm-up rate with tem- 
perature; (2) the heat capacities of any excess liquid 
present, of the metal calorimeter, and of the ‘ Vycor’ glass; 
(3) the heat taken up in evaporation during warming. 

With regard to the first, a typical portion of warm-up 
curve is shown in Figure S. In the diagram, the resistance 
change AR is obtained by a straight line extrapolation 
along DC, C being in the centre of the 30 sec heating 
period. However, a straight line extrapolation was not 
permissible since the rate of warm-up changed rapidly 
with temperature due to the rapidly varying heat capacity 
and the variation in heat leak through the film. Allowance 
had also to be made for the time taken for temperature 
equilibrium to be attained within the calorimeter, the 
initial non-equilibrium being indicated by the curved 
portion C’D. These effects meant that values obtained 
using the simple straight line extrapolation had to be 
corrected by an amount varying from 20 per cent at 1:5°K 
to less than | per cent above 2:1°K. 

The calorimeter contained a small amount of free liquid 
in excess of that inside the ‘ Vycor’ pores. By following the 
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adsorption isotherm during filling, this free liquid could be 
measured accurately, and amounted to 1-70 x 10-° g and 
4-80 x 10-3 g in two experiments, as compared with 0-48 g 
contained in the ‘Vycor’. 

The copper calorimeter weighed 53 g, and the *Vycor’ 
disks when dry weighed together 14-77 g. The heat capac- 
ities were estimated from the specific heats of copper and 
‘Pyrex’ given by Wolcott.!® These amounted together to a 
1 per cent correction at 2:3°K and less at lower tempera- 
tures. 

As the calorimeter temperature T rose and the vapour 
pressure P increased, liquid helium evaporated, thus in- 
creasing the mass of helium gas m, in the dead space in the 
inlet tube. A quantity dQ/d7T was subtracted from the 
apparent heat capacity, using the relation 


dQ L dm dP 

I) Slang 8 as We 
In this equation L was taken as the latent heat of bulk 
liquid helium. In a subsidiary experiment in which the 
“Vycor’ disks were replaced by solid copper disks of the 
same size, 7 was measured as a function of pressure at 
various temperatures enabling dim/dP to be deduced. This 
correction amounted to 10 per cent at 2:3°K falling to 
1 per cent at 1-6°K. 

The final results for the specific heat of the helium within 

the ‘Vycor’, after application of the corrections, are given 
in Figure 6. 
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Figure 6. Specific heat versus temperature of liquid helium-4 contained 
within ‘Vycor’ porous glass. Circles and squares refer to different 
experiments 


Discussion 


It is apparent from Figure 3(a) that the graph of flow 
rate F versus temperature is convex towards the origin, as 
distinct from those in somewhat wider channels or with the 
saturated helium-II film where it is concave. The question 
naturally arises as to where such a curve ultimately reaches 
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the abscissa, if at all. The situation is complicated by the 
effects of the pore size distribution, which presumably 
result in a small amount of flow right up to the lambda 
point of bulk liquid. As a simple criterion we shall assume 
that the region over which (1/F)(dF/dT) assumes its 
maximum value represents the onset temperature. If the 
curve actually met the temperature axis this function would 
be infinite. The graph in Figure 3(b) shows that d(log F)/dT 
reaches a maximum value between 2:0 and 2-1°K for 
specimens X and Y and we shall assume that the onset 
temperature for the predominant pore size occurs within 
this interval. The straight portion of the curve, relating to 
onset of superflow, was found to be similar for different 
samples having onset temperatures down to 1-4°K, being 
simply shifted parallel to the temperature axis with little 
change of slope. 

With regard to the convexity towards the origin of our 
flow rate versus temperature curve, it is worth noting that 
the curves obtained by Brewer and Mendelssohn? ”° for 
the flow of unsaturated films were of very similar shape 
except that at temperatures within 0-1°K of the onset 
temperature they straightened out and intercepted the 
temperature axis at a finite angle. 

This onset temperature for superflow may be compared 
with the temperature 2:06°K at which the specific heat 
curve shows.a maximum, and it appears that the two 
temperatures are identical within about 0-1°K. On the 
other hand, Long and Meyer" have drawn attention to an 
apparent difference of about 0-4°K between the specific 
heat maximum and the onset temperature of unsaturated 
films. However, in arriving at this figure they found it 
necessary to analyse the results of different authors using 
different adsorbents, and to compare readings at the same 
value of P/P,'®11. Adsorption isotherms measured by 
different authors vary considerably even on the same 
material, so the actual thicknesses of film at the chosen 
value of P/P) might have been considerably different. The 
present comparison is free from this difficulty. 

Another point of interest lies in the relationship between 
pore diameter and onset temperature of superfluidity. 
Atkins et al.° have reported results for several specimens of 
‘Vycor’ and these are listed below along with our values in 
Table 2. In order to compare these results with ours it is 
necessary to ascertain that the same criterion for choosing 
the onset temperature has been used in both cases. Atkins 
et al. chose the onset temperature as that at which the liquid 
level moved at about 0-01 cm sec! and using rough figures 
for the dimensions of their apparatus this corresponds to 
about 0:01 cm* min! per cm?. Reference to Figure 3(b) 
shows that this is approximately the same criterion as 
chosen by us. 

It would be of interest next to compare the ‘ Vycor’ pore 
diameter appropriate to a certain onset temperature with 
the thickness of an unsaturated film of the same onset 
temperature. This, however, can only be done approxi- , 
mately because of the uncertainty of film thickness in the 
experiments on unsaturated films. In the fourth column 
of Table 2 we give the percentage of saturation pressure 
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Table 2. Onset Temperatures in ‘Vycor’ and in 
Unsaturated Films 


Onset Average ne ae oe 
e t thickn 
Author ORO | intern (per cent) (Strauss) (eadorihse) 
(A) (A) (A) 
Atkins et al. 1-36 42 48 11 16-5 

Atkins et al. 1-64 50 719 17 28 
Atkins et al. 1-68 57 82 18:5 31 
Atkins et al. 1-78 57 88 23:5 = 
Present work | 2:0-2:1 70 96-98 - = 


necessary to yield the same onset temperature as the 
“Vycor’. These figures were obtained by smoothing the 
results obtained by Brewer and Mendelssohn? and by 
Long and Meyer.’ The fifth and sixth columns list the 
corresponding unsaturated film thicknesses, calculated 
using respectively the isotherms of Strauss! and of 
Frederikse and Gorter.” Of the several different isotherms 
published, these two represent extreme values. We have 
calculated the thicknesses on the assumption that the film 
consists of a helium substrate 4-8 A thick and of density 
0-21 gcm-%, and that this is covered by adsorbate of liquid 
density.12 Unfortunately, the isotherms do not extend to 
high enough pressures for a comparison with the present 
data. 

It thus appears that the film thickness at which, at a 
given temperature, superfluidity sets in is considerably less 
than the corresponding average pore diameter. If the 
porous structure consisted of large cavities connected by 
narrow constrictions, then it is possible that the actual 
diameter determining the flow in the ‘ Vycor’ would be less 
than the mean pore diameter. However, the specific heat 
isa bulk property not affected appreciably by constrictions, 
and if the structure were as indicated above one would 
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expect the onset temperature of superflow to be below that 
of the specific heat maximum. Our results thus favour the 
alternative suggestion that depression of the lambda 
transition is caused by the forces of attraction to the wall. 
Since in ‘ Vycor’ glass the liquid is completely surrounded 
by solid walls, the effect of attractive forces is greater than 
in the film, where the forces operate on only one side. 


We are grateful to Dr. M. E. Nordberg of the Corning 
Glass Research Laboratories for supplying us with the 
glass samples used in this research. 
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Discussion on Enthalpies of Mixtures 


M. Rahemann  Petrocarbon Developments Lid, Manchester 


Received 17 August 1960 


ComMISSION I of the International Institute of Refrigera- 
tion, meeting at Eindhoven at the end of June, devoted a 
special session under the chairmanship of the author to 
‘Enthalpy Ignorance’, the lack of information on the 
enthalpies of gaseous and liquid mixtures. Short contri- 
butions were made by invited speakers followed by a 
general discussion of the subject. 

In an introductory statement, Dr. P. J. Haringhuizen 
(Staatsmijnen, Limburg) pointed out that the design of low 
temperature gas separation plants requires a knowledge of 
the enthalpies of gaseous mixtures as a function of pres- 
sure, temperature, and composition. Whereas fairly re- 
liable charts and tables exist for the thermodynamic 
properties of most of the components occurring in these 
mixtures, the properties of the mixtures themselves are 
little known. 

Atlow pressures and well above the critical temperatures 
of all the components, the enthalpy may, as a first approxi- 
mation, be taken as a linear combination of the enthalpies 
of the constituents, neglecting the heat of mixing. However, 
at high pressures and low temperatures, this is certainly 
not permissible. As an example, Dr. Haringhuizen showed 
a plot of the enthalpy of a 3:1 hydrogen-nitrogen mixture 
as a function of pressure up to 400 atm for several tem- 
peratures above O°C, which had been established 
experimentally by Michels and colleagues at Amsterdam. 
On the same plot calculated enthalpies were shown, using 
a linear combination of the component enthalpies based 
on the total pressure and, alternatively, on the partial 
pressures of the components in the mixture. Both theo- 
retical curves showed marked deviations from the experi- 
mental values, the curves based on the total pressure being 
‘worse’ than those using partial pressures, and the 
deviations became greater the lower the temperature. 

An additional difficulty arises when, as is frequently the 
case at the required temperatures and pressures, one or 
several of the pure components exist by themselves in the 
liquid state only. It is evidently incorrect to insert the 
enthalpies of the pure liquids in the linear combinations 
intended to give the enthalpies of the gaseous mixtures. 

Proposals have been made to calculate enthalpies of 
hydrocarbon mixtures on the assumption that they depend 
only on the molal average boiling point, but these methods 
are largely empirical and have no sound thermodynamic 
basis. 
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In addition to the enthalpies of gas mixtures, those of 
the corresponding liquids are also required for design 
purposes, and our knowledge of the heats of mixing of 
liquids is still very fragmentary. Moreover, little or no 
information is available to assess the heat to be withdrawn 
from a mixture under conditions of partial condensation 
in the presence of non-condensing components. 

Dr. Z. Dokoupil (Kamerlingh Onnes Laboratory, 
Leiden) gave an example of how the enthalpy of a gas 
mixture can in principle be evaluated on the basis of 
P-V-T data for the individual components with the help 
of a suitable equation of state. Taking the 3H,+ N, mix- 
ture treated by Dr. Haringhuizen, for which experimental 
data are available, and using the Beattie-Bridgeman 
equation of state, the relevant constants of which have 
been tabulated for the pure components, he calculated the 
constants for the mixture with the conventional combina- 
tion rules. He then computed the enthalpy of the mixture 
along the 100°C isothermal, using well known thermo- 
dynamic relationships, and found that the results hardly 
deviated from the curve based on the experimental P-V-T 
data up to 400 atm pressure. If the constants of the 
Beattie-Bridgeman equation are adjusted to fit the 
experimental data for the pure components as accurately 
as possible in the field of temperature and pressure 
covered, still better agreement can be obtained for the 
mixture. 

It may be inferred that this method is applicable to 
certain gas mixtures at relatively high temperatures up to 
pressures of several hundred atmospheres, but it may well 
become inaccurate at low temperatures and cannot be 
used below the critical temperatures of any of the com- 
ponents. 

Dr. F. Din (British Oxygen Research and Development 
Ltd., London) approached the problem from the opposite 
side by demonstrating a method of calculating the heat of 
mixing of liquid oxygen and nitrogen at low pressure from 
liquid—vapour equilibrium data. This system had been 
investigated experimentally by the author over the com-. 
plete range of compositions and an adequate field of 
temperature and pressure. The results enabled the expres- 
sion In (y,/yz) to be calculated from the relative volatility 
and the fugacities of liquid and vapour, where y, and y, 
are the activity coefficients of nitrogen and oxygen. It will 
be recalled that the activity coefficient may be defined by 
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Si =i yi Xi, where f, is the fugacity of component i in the 
solution, f; the fugacity of the pure component at the same 
temperature and at the total pressure of the solution, and 
x; is the mole fraction. 

Now the excess free energy AG“) (free energy of real 
solution minus free energy of ideal solution) can be 
evaluated from the equation 


dAG®) 
(2°), 
x; ip Y2 


and this was done for temperatures of 80 and 90° K. 

The required heat of mixing is in fact the excess enthalpy 
AH™®) of the mixture since the enthalpy of mixing of the 
ideal solution is zero. 

Unfortunately, the corresponding integral equation 
involving the excess enthalpy of mixing AH‘®) does not 
allow the latter to be evaluated unless it is a very simple 
function of composition, e.g. of the form Ax, x2, where A 
is a constant. Values obtained with the Gibbs—Helmholtz 
equation are not sufficiently accurate. It may, however, be 
tentatively assumed that liquid mixtures of oxygen and 
nitrogen behave as ‘regular’ solutions for which, by 
definition, the excess entropy of mixing AS‘) = 0. In this 
case, since 

AG®) = AH® =TFTAS® 
then 
AH® ="AG® 


Assuming that AH has the form Ax, Xo, it was possible 
to determine A from the integral equation mentioned 
above and, hence, AH‘) for | atm. The resulting curve is 
very close to the curves for AG‘) at 80 and 90° K and it can 
therefore be assumed that the assumption of a regular 
solution is justified. 

Professor W. C. Edminster (Oklahoma State University) 
reviewed the work carried out in recent years and planned 
for the immediate future in the U.S.A. to overcome the 
‘enthalpy ignorance’, the significance of which had been 
fully realized in that country where, it was estimated, it 
was costing over $1,500,000 annually. 

Enthalpies of gaseous and liquid mixtures can be com- 
puted by way of thermodynamic relations from accurate 
P-V-T-X data or determined experimentally by calori- 
metric measurements. The results can be used to develop 
compilations of partial enthalpies as functions of tem- 
perature, pressure, and composition for each component. 
The difficulty is to find a suitable composition parameter. 

The former method was used twenty years ago by 
Benedict and colleagues, who developed an 8-constant 
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empirical equation of state, correlated the resulting 
fugacity coefficients for light hydrocarbons, and used these 
for calculating enthalpies of mixtures. The same fugacuy 
correlations were later used by Edminster and Canjar to 
compute partial enthalpies. An attempt was made to 
generalize these in terms of reduced parameters, but 
extrapolation to obtain partial enthalpies of other sub- 
stances in different mixtures was found to be impossible. 

The correlation of AH values computed for specific 
systems directly from P-V-T-X data appears more 
advantageous than the reverse procedure. 

A project to obtain enthalpies of mixtures has been 
initiated in the U.S.A. by the Natural Gasoline Associ- 
ation. This project, which will take about 8 years to 
complete and will cost $300,000, will involve experimental 
work at two or three universities, data correlation, theo- 
retical work, and preparation of charts, tables, equations, 
and computer programmes. 

Experimental work has been started at a number of 
places, including the Chemical Engineering Department 
of the University of Michigan, the Institute of Gas Tech- 
nology at Chicago, and in the Chemical Engineering 
Laboratory at Oklahoma State University. Here Professor 
Edminster is intending to use an isobaric flow calorimeter, 
an isothermal throttling calorimeter, and an isobaric 
calorimeter to measure integral heats of vaporization. 

Professor G. G. Haselden (Leeds University). consider- 
ing ways of determining enthalpies of mixtures, stressed 
the importance of direct methods in view of the magnitude 
of the task. The measurement of Joule-Thomson co- 
efficients provided a quick.and reasonably direct method 
for single-phase systems. Heats of mixing and flow calor- 
imeter methods were suitable in the two-phase region. 
The boiling type calorimeter tended to be inaccurate for 
enthalpy determinations, especially at high pressures. 
owing to the instability of the boiling process and the large 
thermal capacity of the apparatus. Work on heats of 
mixing had been started at Leeds University and a con- 
densing type flow calorimeter had been built at Imperial 
College, London. 

Summing up the contributions and the subsequent 
discussion, the Chairman proposed that Commission I of 
the International Institute of Refrigeration should take 
steps to-stimulate work on the enthalpies of mixtures, both 
experimental and theoretical, the latter to include funda- 
mental work as well as the correlation of results in a form 
suitable for design purposes. Co-ordination of work in 
progress should be encouraged both within individual 
countries and on an international basis. A motion embody- 
ing these proposals was carried by the meeting. 
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Letters to the Editors 


A Proposal for a New Nuclear Demagnetization Method 


KurtTI, Simon, Robinson, and Spohr succeeded in 1956 
in the difficult experiment of cooling nuclei of metallic 
copper to a temperature of 10-°°K. Later on, their method 
was developed further and new results were published.! 
The method they used had been proposed some years 
previously by Gorter? and by Kurti and Simon.* 4 It 
consists of a two-stage demagnetization. 

First an electronic paramagnetic salt (pill A) is cooled 
to about 0-01°K by ordinary adiabatic demagnetization. 
This heat sink is connected thermally with a pill B con- 
taining nuclei with a magnetic moment. The pill B is placed 
in a strong magnetic field, but the field at the electronic 
stage must be weak. It follows that there must be a fairly 
long connection between the pills. After removal of the 
heat of magnetization from pill B to pill A, the magnetic 
field is reduced and the substance of pill B cools down. 

In order to overcome the difficulty of the heat transport 
between the pills, Kurti et al. chose copper for the material 
of the link and also for pill B. 


Trig. 
axis 


BS Figure 1. Principle of the arrangement 
with core C and shell I 


It would be advantageous if protons could be taken as 
the nuclear stage instead of copper nuclei. Protons have 
a magnetic moment which is 3-8 times that of both the 
copper isotopes. Thus the entropy reduction would be 
15 times larger at the same initial field and temperature, 
as this reduction is proportional to the square of the mag- 
netic moment. However, the difficulty then is the heat 
transfer from the nuclear to the electronic stage. 

The problem might be solved by the application of 
monocrystals with an anisotropic g factor, such as a 
number of paramagnetic salts of the rare earths group 
possess. With these crystals a low temperature can be 
achieved by rotation of the magnetic field.® This property 
suggests the idea of a pill built up of a core containing the 
nuclear stage (protons) surrounded by a shell made of one 
of the above-mentioned paramagnetic rare earth crystals 
(Figure 1). Now the nuclear substance is on all sides in 
heat contact with the electronic substance. 

Cerium magnesium nitrate is a suitable substance for 
shell I. It has a trigonal structure, and the g factor perpen- 
dicular to the trigonal axis (g,) is much larger than the g 
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factor parallel to this axis (g,). Very low temperatures 

(0-01°K) can be obtained by magnetizing this salt in a field 

perpendicular to the trigonal axis, and demagnetizing it 

by rotation of the field to a direction parallel to the trigonat 
axis. 

The experimental procedure, in principle, is then as 
follows: 

(1) The pill is magnetized perpendicularly to the trigonal 
axis with the strongest possible field at the lowest 
temperature attainable by pumping liquid helium. 

(2) The heat contact between pill and bath is broken. 

(3) The magnetic field is turned slowly in the direction of 
the trigonal axis, resulting in a reduction of the tem- 
perature of the shell I, which cools in its turn the core C. 

(4) The magnetic field is reduced to zero and the core cools 
down. 


Trig. 
axes 


C: Lanthanum magnesium nitrate I: Cerium magnesium nitrate 
Il: Neodymium magnesium nitrate 


Figure 2. Two different examples consisting of core C, shell I, guard 
shell IT, and coil A 


We intend to use a pill with a guard ring. Two different 
examples are given in Figure 2, showing cross-sections 
parallel to the (vertical) cryostat axis. The core C is made 
of a monocrystal of, for example, lanthanum magnesium 
nitrate. The protons of the crystallization water will be 
magnetized. Shell I, which must have as good a heat 
contact as possible with core C, is made of one or more 
monocrystals of cerium magnesium nitrate. A fair heat 
contact can be achieved by grinding the outer surface of C 
and the inner surface of the shell I, and by applying some 
silicon grease or other cement between them.*-8 

Shell II is made of one or more monocrystals of neo- 
dymium magnesium nitrate with the trigonal axis parallel 
to that of the cerium salt and mounted in such a way as to 
have bad heat contact with shell I. As g, <g, for this salt 
also, it too is cooled by rotation of the field. The specific 
heat of the neodymium salt is much higher than that of the 
cerium, so it can serve well as a guard shell. 

In order to measure the temperature of the core C, a 
coil A is provided on C and thus inside shell T; it is the 
secondary coil of a Hartshorn bridge. 
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The heat transfer of the shell I to the lattice of core C 
can be made fairly high. That of the lattice to the protons 
of the core is very low in accordance with the very long 
lattice—nuclei relaxation time. This relaxation time, how- 
ever, can be reduced by the incorporation of paramagnetic 
ions in the core C as indicated by Bloembergen,? by Rollin 
and Hatton,?° and by Khutsishvili!! (either Cr3+ in the 
alum, or neodymium in the nitrate). 

These paramagnetic contaminations will have a large 
effect on the susceptibility of the core. Kurti,!2 however, 
showed that below about 0:02°K the susceptibility of the 
impurities becomes independent of temperature, so that 
below that temperature any change in susceptibility can be 
ascribed to the nuclear spin system. 

_ It is hoped that Kurti’s conclusion will also remain 
substantially correct in our experiments, although the 
amount of impurities will be much higher. 

If the longitudinal relaxation time ¢, of protons or other 
nuclei in an anisotropic (eventually diluted) salt such as 
cerium or neodymium magnesium nitrate were of the 
order of several minutes at temperatures of about 0-01°K, 
a much simpler procedure for the nuclear cooling could be 
followed. Then, after slow rotation cooling of the salt, 
one should have to wait only a reasonable time for the 
electron system to cool the proton system. Following this, 
the demagnetization would cool the nuclei to a very low 
temperature. The warming up of the nuclear by the elec- 
tronic system would be slow enough to allow the tem- 
perature measurement. 

‘A more detailed account of this proposal will appear in 
the proceedings of Commission I of the International 
Institute of Refrigeration, held at Eindhoven on 28-30 
June 1960. 


B. S. BLAISSE 
Laboratorium voor Technische Fysica, 
Delft, 
Holland. (18 July 1960) 
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The Cooling of Specimens for Calorimetry at Very Low 
Temperatures 


ONE very tricky problem in calorimetry at extremely low 
temperatures is the cooling of a specimen which, for the 
measurement, is thermally insulated inside a vessel such 
as V in Figure 1, the vessel itself being submerged in, for 
example, a bath of liquid helium. This cooling may be 
effected either by temporarily introducing an exchange 
gas, such as helium, into V, or by making a thermal con- 
tact between the specimen E and the wall through a 
mechanical process. 

The exchange gas method has been generally used for a 
long time and yields excellent results with metals and 
alloys.1 However, it has the disadvantage of requiring 
pumping for a considerable time after the right tempera- 
ture has been reached. The vacuum channel must have a 
very small diameter in order to prevent radiation heat 
from entering and, moreover, the exchange gas is ad- 
sorbed on the specimen. If the specimen has a large area 
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(powdered or sintered products), the quantity of gas so 
adsorbed may be considerable and may cause distur- 
bances which are not well understood. If measurements are 
made at very low temperatures, where the specific heat of 
the solid is itself considerably reduced, the calorific capacity 
of the adsorbed fluid, which varies only little as a function 
of the temperature, becomes an important correction. 

In order to be efficient, the mechanical contact requires 
considerable pressure, as Berman? has shown in his work 
on the effect of contact force on conductance. Sometimes 
the force is transmitted directly from the outside, but 
preferably it is amplified at the low temperature. The 
most common mechanical method consists of locking 
the specimen in a clamp kept at the temperature of the 
vessel V controlled by vertical traction.* Because of the 
clamp, and in spite of the amplification, the size of the 
control wire used remains fairly important, bringing an 
amount of heat which is not negligible. Also, since the 
clamp is necessarily rather massive, a fairly long time 
elapses before its temperature is reduced to that of the 
vessel V. 

Finally, as it is required that, at the moment when the 
contact is broken, no vibrations which might generate 
heat are created in the specimen, it is necessary to break 
the contact very gradually, the solid frictions in the trans- 
mission are therefore undesirable. 

In our method, described in the following, a hydraulic 
contact control is used: the force is exerted by means of a 
metal bellows S, made of tombac, which is simply inflated 
with helium. The filling is effected through a capillary 
t-t’. The force which can be applied is of the same order 
of magnitude as that applied with a clamp. With the 
tombac bellows used by us, with an inner diameter of 25 
mm, and at a pressure of 4 kg/cm?, the force is 20 kg. 
Under the effects of the pressure, the bellows expands and 
pushes the specimen E, suspended on the nylon threads 
F and: F’, against the plate P at the upper part of the 
frame C. In the case shown in the drawing, C has the 
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same temperature as V. The descent of the bellows is 
controlled by a fine regulating valve, r, and the move- 
ment of a few millimetres which will break the contact 
can be effected as slowly and progressively as may be 
desired. 

Without taking any special precautions regarding con- 
tacting surface at P and P’ or on the specimen, we have 
already been able to effect breaks entailing a liberation of 
heat of only 100-200 ergs using this method. That is, the 
heat release is of the same order of magnitude as that 
obtained with the most favourable contact surfaces by the 
mechanical methods used previously (several hundred 
ergs). As the capillary is located in liquid helium, all 
thermal disturbances are eliminated. Finally, the device 
is far simpler than the mechanical arrangements proposed 
hitherto. 

This cooling method was of particular use in some work 
we are carrying out on sintered beryllium oxide which, 
being an insulator with a particularly high Debye 
temperature, has a very low specific heat: a 65 g specimen 
at 2° K has a calorific capacity of only 4 x 10~® joules per 
degree. 


J. ASLANIAN 
L. WEIL 


(26 September 1960) 


Université de Grenoble, 
France. 
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Book Reviews 


Cryophysics. K. Mendelssohn. ‘Tracts on Physics and 
Astronomy’ Series. (Interscience Publishers, 1960) 192 
pages. Cloth $4.50 (34s.); Paper $2.50 (19s). 


Within the past five years the literature of low temperature 
physics has grown at a remarkable rate. In addition to 
review articles on particular topics in the new Handbuch 
der Physik and in Advances in Solid State Physics, we now 
have two ‘Progress’ séries devoted specifically to the field, 
each of which has so far produced two volumes, and 
several books dealing with liquefiers and low temperature 
techniques. To supplement these, voluminous conference 
reports appear at regular intervals. But for the under- 
graduate or the beginner in research who seeks a broad 
general survey of the subject, none of these is suitable, at 
any rate for first reading. The latest edition of Jackson’s 
excellent monograph appeared in 1955; much has hap- 
pened since then, and there has been a real need for a more 
up-to-date introductory text. Dr. Mendelssohn’s new book 
has been written in an attempt to fill this need. Within the 
space of 180 pages, he has managed to touch on practically 
every topic in the field, including many of very recent 
development. The chapter headings follow roughly the 
traditional pattern—cryogenics, thermometry, specific 
heats, magnetism, transport phenomena, superconduc- 
tivity, and helium—and sections are included on nuclear 
orientation, dielectric heat conduction, cryotrons, super- 
conductivity theory (with a rather lukewarm reference to 
the Bardeen—Cooper-—Schrieffer theory), turbulence in 
liquid helium, and liquid helium-3; an additional chapter 
entitled ‘ Miscellaneous’ covers the melting curve, thermal 
expansion, mechanical properties, trapped free radicals, 
bubble chambers, and masers. Each chapter ends with a 
useful list of references to recent review articles. 

It has not been possible to cover this wide range of topics 
without some sacrifice, of course. As the author states in 
his preface ‘detailed treatment has been avoided through- 
out’, but to achieve even a superficial survey in the space 
available he has had to write with the greatest economy. 
At many points a few extra words would have added 
greatly to the clarity, and in some places the beginner is 
likely to be seriously puzzled. For example, on page 21, in 
discussing the Oxford work on nuclear demagnetization 
it is remarked in passing that ‘the rate at which the meas- 
ured nuclear susceptibility decreased again was high’, but 
the implications of this remark remain obscure until 
page 26. Again, on page 62, the sudden appearance of the 
word ‘ferrimagnetism’, once and without explanation, is 
likely to confuse rather than enlighten the student, unless 
he supposes it to be a misprint. On page 44, in the section 
on electronic specific heat, it is not made clear that equa- 
tion (3-4) for an ideal Fermi—Dirac gas applies only when 
the departure from classical behaviour is very small, and 
is quite irrelevant to the subsequent discussion of electrons 
in a metal. And on page 168, the student is liable to be 
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puzzled by the dimensionless ‘energy gaps’ A and 6. 
These obscurities, and others, are possibly inevitable when 
so much is compressed into the space available, but they 
make it difficult to recommend the book unreservedly to 
the beginner. For those who already have reasonable 
familiarity with the subject, however, it can be recom- 
mended as a useful summary of recent developments. 


R. G. CHAMBERS 


Liquid Helium. K. R. Atkins. ‘Monographs on Physics’ 
Series, (Cambridge University Press) 312 pages. 60s. 


Most of the distinctive characteristics of liquid helium-4 
were essentially known fifteen years ago or more, and 
much of the basic theory for understanding them laid 
down. Strong elements of mystery and uncertainty were 
subsequently attacked with vigour, but it was only a few 
years ago that a theoretical picture emerged, supported 
by experiment. which promised to lead to a completely 
satisfactory solution. Professor Atkins’ book is thus well 
timed, coming as it does shortly after this great advance 
in the subject. While there are undoubtedly many interest- 
ing discoveries still to be made (and some have been made 
since the book was written), the approach seems now to 
be rather one of a mopping-up operation than a funda- 
mental attack. As the author believes, the ultimate solution 
will probably be along similar lines to those he describes, 
and it will thus be some time before his book becomes 
seriously out of date. 

The book inclines to explain the current situation 
rather than give an historical exposition, an approach 
which has its advantages in a subject with a long and 
sometimes confusing history. It makes easy and stimu- 
lating reading, and succeeds in giving a wide survey of 
work on helium-II, although naturally one may not 
always agree with the emphasis placed on the various 
aspects. The treatment of thermodynamics and vortex 
line theory forms a particularly useful introduction to the 
subject, although the impression may sometimes be re- 
ceived that the observations are better understood than 
is in fact the case. Two chapters on liquid helium-3 
and on liquid helium-3—helium-4 mixtures are also in- 
cluded. Fundamental discoveries about helium-3 are 
being made at a rapid pace and this section has 
already become somewhat out of date. It still provides a 
most useful and interesting introduction because of the 
author’s clear exposition. The book as a whole will 
certainly become compulsory reading for research 
students—it is the only book on the subject of which 
this can be said—and a stimulating, if not quite 
exhaustive, reference for others. . 

1D. F. BREWER 
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